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Summary. — Optics and interferometry with matter waves is the art of coher¬ 
ently manipulating the translational motion of particles like neutrons, atoms and 
molecules. Coherent atom optics is an extension of techniques that were developed 
for manipulating internal quantum states. Applying these ideas to translational 
motion required the development of techniques to localize atoms and transfer popu¬ 
lation coherently between distant localities. In this view position and momentum are 
(continuouse) quantum mechanical degree of freedom analogous to discrete internal 
quantum states. In our contribution we start with an introduction into matter-wave 
optics in section 1, discuss coherent atom optics and atom interferometry techniques 
for molecular beams in section 2 and for trapped atoms in section 3. In section 4 we 
then describe tools and experiments that allow us to probe the evolution of quantum 
states of many-body systems by atom interference. 
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1. Optics and interferometry with atoms: an introduction 

Interference is one of the hallmark features of all wave theories. Atom interferometry 
[1] is the art of coherently manipulating the translational and internal states of atoms 
(and molecules), and with it one of the key experimental techniques to exploit matter 
waves. From the first experiments demonstrating the wave-like nature of light [2, 3] 
to the ground-breaking achievements of matter-wave interferometry with electrons [4], 
neutrons [5], atoms [6] and even large molecules [7], interference has led to new insights 
into the laws of nature and served as a sensitive tool for metrology. 

Interference with atomic and molecular matter-waves now forms a rich branch of 
atomic physics and quantum optics. It started with atom diffraction of He from a LiF 
crystal surface [8] and the separated oscillatory fields technique [9] now used in atomic 
clocks. Broadly speaking, at the start of the 20th century atomic beams were developed 
to isolate atoms from their environment; this is a requirement for maintaining quantum 
coherence of any sort. In 1924 Hanle studied coherent superpositions of atomic internal 
states that lasted for tens of ns in atomic vapors [10]. But with atomic beams, Stern- 
Gerlach magnets were used to select and preserve atoms in specific quantum states for 
several ms. A big step forward was the ability to change internal quantum states using 
RF resonance, as demonstrated by Rabi et al. in 1938 [11]. Subsequently, long-lived 
coherent superpositions of internal quantum states were created and detected by Ramsey 
in 1949 [9] which is the basis of modern atomic clocks and most of quantum metrology. 

Applying these ideas to spatial degrees fo freedom required the development of tech¬ 
niques to transfer atoms coherently between different locations. The simplest way is to 
create coherent superpositions of states with differetn momenta. Here, coherently means 
with respect to the phase of the de Broglie wave that represents this motion. 

We give an introduction ot the coherent atom optics techniques in sections 1 and 2 
(A much more complete overview is given in the review by Cronin et al. [1]). Although 
some experiments with Bose-Einstein condensates are included, the focus of these two 
first sections is on linear matter wave optics where each single atom interferes with itself. 
Techniques for trapped atoms are then discussed in section 3 and in section 4 we describe 
recent tools and experiments, where atomic interference is used to probe the complex 
quantum states of interacting many-body systems. 

IT. Basics of matter-wave optics. - In this first section we discuss the basics of 
matter-wave optics and illustrate the similarities and differences to the more familiar 
light optics. This is by no way a detailed and in depth theoretical discussion, but should 
merely highlight the differences and similarities between matter-waves and light. For a 
detailed theoretical discussion we refer the reader to the lectures of Ch. Borde in these 
proceedings. 

IT.l. The wave equations. A first approach to comparing matter-wave optics to light 
optics is to study the underlying wave equations. The differences and similarities between 
light optics and matter-wave optics can then be nicely illustrated in the following way: 
Light optics is described by Maxwell’s equations. They can be transformed and rewritten 
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as the d’Alembert equation for the vector potential A 


( 1 ) 



(? dt"^ 


A(r, t) = 0. 


Here for simplicity, we have assumed an isotropic and homogeneous propagation medium 
with refractive index n, and c denotes the speed of light in vacuum. 

If we now consider a monochromatic wave oscillating at an angular frequency u, i.e. 
we go to Fourier space with respect to time, A(r, t) —)■ A(r, w) exp(—zwt), we obtain the 
Helmholtz equation for A(r,a;) 


( 2 ) 



A(r, w) = 0. 


The propagation of matter-waves for a non-relativistic particle is governed by the 
time-dependent Schrodinger equation. For non-interacting particles, or sufficiently dilute 
beams, a single particle approach is sufficient: 


( 3 ) 




,, , .d'lblY.t) 

= ih ———, 


where m is the mass of the atom and V (r) is a scalar potential. 

Comparing equations (1) and (3) one finds a hyperbolic differential equation for light 
optics (equation (1)), whereas the Schrodinger equation (3) is of parabolic form. This dif¬ 
ference in the fundamental wave equations would suggest significantly different behavior 
for matter-waves and light waves. 

If one only considers time independent problems, like, for example, propagation of a 
plane wave in a time independent potential V{r), we can eliminate the explicit time de¬ 
pendence in equation (3) by substituting —>■ ■i/'(r) exp(—where E is the total 

energy, which is a constant of the motion for time independent interactions. The prop¬ 
agation of the de Broglie waves is then described by the time independent Schrodinger 
equation 


( 4 ) 




tp{r) = 0. 


We can define the local A:-vector for a particle with mass to in a potential V (r) as 
(5) k{r) = ^^/2m{E - C(r)), 

and equation (4) becomes equivalent to the Helmholtz equation (2) for the propaga¬ 
tion of electromagnetic fields. Therefore, at the level of equations (2) and (4), for a 
monochromatic and time-independent wave, matter waves and classical electromagnetic 
waves behave similarly. If identical boundary conditions can be realized the solutions for 
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the wave function in matter-wave optics and the electric field in light optics will be the 
same. Many of the familiar phenomena of light optics, like refraction (see section 1'4) 
and diffraction (see section 2’1), will also appear in matter-wave optics. 

1’1.2. Dispersion Relations. The dispersion relation of a wave relates its energy to its 
/c-vector. Dispersion relations become apparent in the Helmholtz equations (2) and (4) by 
applying a Fourier transform to the spacial coordinates of the fields, i.e. by substituting 
'tp{r) — )■ ^/;(k)exp(tk • r). Writing ko as the wave-vector in vacuum, we obtain vacuum 
dispersion relations which are linear for light 

(6) oj = — = cfco, 

n 

and quadratic for matter-wavesi}) 

_h^_ hk^ _ E-V 
^ 2m 2mn^ h 

An important fact we observe in equation (7) is that the mass enters the dispersion 
relation for matter-waves. Moreover, the quadratic dispersion relation for matter-waves 
causes even the vacuum to be dispersive. A consequence of this dispersion is, for example, 
the spreading of a wave packet, which happens even in the longitudinal direction. For 
example, a Gaussian minimum uncertainty wave packet, prepared a time < = 0 as 

(8) V’(i-,t) ocy A(k,ko)e*(’^ '-“‘)dk, 


with a momentum distribution 


( 9 ) 


(k-kp)^ 

A(k, ko)oce 


spreads even when propagating in vacuum, that is in the absence of a refractive medium. 
Here, ax{0) = h/uk is the width of the wave packet in position space. For our Gaussian 
minimum uncertainty wave packet the spreading results in a time dependent width in 
real space which can be written as 


(10) ax{t) = ^cr2(0)-f 0-2(0) 

This spreading of the wave packet is nothing else than the wave mechanical equivalence 
of the dependence of the propagation velocity on the kinetic energy of a massive particle 


(^) In a relativistic description the dispersion relation is given by —h which 

reduces to a; ~ —h 2 mn^ non-relativistic limit . The difference with equation (7) is 

caused by the energy associated with the rest mass of the particle. 
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in classical mechanics. In wave mechanics, the wave packet spreads due to its different 
/c-space components moving at different velocities. 

I’l.S. Phase and Group Velocity. Another consequence of the different dispersion 
relations is that the phase velocity Vph and group velocity Vg for de Broglie waves are 
different from those of light. In a medium with refractive index n (section 1'4) one finds 
for matter-waves 


( 11 ) 

'^ph 

( 12 ) 





V 

2 ’ 


1 

-Vq = V. 

n 


The group velocity Vg, as given in (12), corresponds to the classical velocity v of the 
particle. For a wave packet this corresponds to the velocity of the wave packet envelope. 
For matter-waves the vacuum is dispersive, that is Vph 7 ^ Vg for n = 1. Furthermore, Vph 
and Vg are both inversely proportional to the refractive index(^), with 

(13) VpHVg = ^=vy2. 

It is interesting to note that similar phenomena, like non-linear dispersion relations, 
the spreading of wave packets and the non-equivalence of phase and group velocity can 
also be found in the propagation of electromagnetic waves in refractive media, or in wave 
guides. The details of this correspondence depend on the detailed dispersion character¬ 
istics in the refractive index n{k) of the medium, or the wave guide. 

1’2. Path Integral Formulation. - The wave function description of the propagation 
of light or matter waves is very illustrating and powerful. Nevertheless many problems 
can be solved more easily by an equivalent approach developed by Feynman, where 
the amplitude and phase of the propagating wave at a position in space and time are 
expressed as the sum over all possible paths between the source and the observation 
point(^). This method can, for most cases in matter-wave optics, be simplihed and gives 


(^) In a relativistic description the phase and group velocity for matter-waves can be obtained 
as follows {E is now the total energy including the rest mass): 

Lj IE me? hko <? u 
k n hko nhko 2nm v 2 ’ 

doj hke^ 1 hkoc^ 1 hko 1 

dk E n E n m n 

where Vph and Vg are again both inversely proportional to the refractive index, but the product 
VphVg is now given by VphVg = (? , which is the same as for light. 

(®) A good and easily readable summary, adapted for atom optics is given by P. Storey and 
C. Cohen-Tannoudji in ref. [12]. 


'^ph — 

^9 = 
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straightforward, easily interpretable results. 

In the regime of classical dynamics the path a particle takes is determined by the 
equation of motion. The actual path taken can be found by the principle of least action 
from the Lagrangian 

(14) L=^mz‘^{t)-V{z)=pz-H. 


Here, z is the spatial coordinate, H the Hamiltonian, and the momentum p is defined as 
P ~ classical action is defined as the integral of the Lagrangian over the path T 


(15) Sr = [ L{z{t),z(t))dt. 

Jti 


In general, the dynamics of the system is described by the Lagrangian equations of motion 


(16) 


dL d dL 
dz dt dz 


which are the differential form of the principle of least action and completely equivalent 
to Newtons equations. Note that the principle of least action which defines the classical 
paths for particles is equivalent to Fermat’s principle for rays in classical light optics. 

For a quantum description one has to calculate the phase and amplitude of the wave 
function. As it was pointed out by Feynman [13, 14], the wave function at point a can 
be calculated by superposing all possible paths that lead to a. In general the state of 
the quantum system at time tb is connected to its state at an earlier time ta by the 
time-evolution operator U(ti,,ta) = exp{—H{tb — ta)/h) via 


(17) IV’(4)) = U {tb,ta)\'>P{ta)) ■ 

The wave function 'jp{xbjtb) at point Xb is given by the projection onto position 


(18) 


'tp{xb,tb) = {Xb\^l^{tb)) = 


dXaK{xb, h; Xa, ta)lp{Xa, ta), 


where the quantum propagator K{xb,tb',Xa,ta) is defined as 
(19) K {xb,tb, X a,ta) = {a^h\U {tb,ta)\^a) • 

Equation (18) is a direct manifestation of the quantum mechanical superposition princi¬ 
ple and shows the similarities of quantum mechanical wave propagation to the Fresnel- 
Huygens principle in optics: The value of the wave function at point {xb,tb) is a super¬ 
position of all wavelets emitted by all point sources {xa,ta)- 

Furthermore, the quantum propagator has some properties which are very useful for 
real calculations. One such property comes from the fact that the evolution of a quantum 
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system from time ta to time tb can always be broken up into two pieces at a time tc with 
ta < tc < tb- The calculation can be done in two steps from ta to time tc and then 
from tc to tb using the identity U(tb,ta) = U(tb,tc)U{tc,ta)- Therefore, the composition 
property of the quantum propagator is given by 

( 20 ) K{^Xb-,tb^ Xajta') — {Xb\U{tb-,tc)U{tci^a^\^a^ 

— dXc K(^Xb^ tb^ Xc^ tc^^i^Xc^ tc^ Xa-, ta^- 

This shows that the propagation may be interpreted as summation over all possible 
intermediate states. It is also interesting to note that this composition property applies 
to the amplitudes and not the probabilities. This is a distinct feature of the quantum 
evolution, which is equivalent to the superposition of the electric fields in optics. 

Based on this composition property of the quantum propagator we can give Feynman’s 
formulation of K{xb,tb;Xa,ta) as a sum over all contributions from all possible paths 
connecting [xa^ta) to {xb,tb) [13, 14] 

( 21 ) Kixb,tb;Xa,ta)=NY,e"^^^'", 

r 

where is a normalization and is the sum (integral) over all possible paths connect¬ 
ing (xa,ta) to (xb,tb). Each path contributes with the same modulus but with a phase 
factor determined by Sr/h where S'r is the classical action along the path T. Feynman’s 
formulation is completely equivalent to the formulation of equation (19). 

In the quasi-classical limit, where Sr 3> h, the phase varies very rapidly along the 
path and most of the interference will be destructive, except where the classical action 
has an extremum. Only paths close to the classical path described by equation 16 will 
then contribute significantly to the sum in equation (21). 

The method of path integrals is a very powerful method to solve the problem of 
propagating matter waves. However it is, even for very simple geometries, very hard 
to implement in its most general form. In most cases of matter-wave optics we can use 
approximations to the full Feynman path integral formulation. The possible approxima¬ 
tions follow from the observation that the largest contribution to the path integral comes 
from the paths close to a path with an extremum in the classical action S'r. 

1’2.1. JWKB approximation. The first approximation that can be done is the JWKB 
approximation(^), often also called the quasi-classical approximation: one uses the classi¬ 
cal path to calculate phase and amplitude of the wave function at a specific location, i.e. 

(22) V’ = ^ 

Vp Vp 


(^) This method was first introduced by Lord Rayleigh for the solution of wave propagation 
problems. It was then applied to quantum mechanics by H. Jeffreys (1923) and further developed 
by G. Wentzel, H.A. Kramers and L. Brillouin (1926). 
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In the JWKB approximation one easily sees how the wavefronts and the classical trajec¬ 
tories correspond to each other. For a fixed energy E, a wavefront is given by the relation 
S'(x) = 5'o = constant. One can show(^) that for a scalar potential the wavefronts are 
orthogonal to the classical trajectories. 

In the case of a vector potential A one has to replace the classical momentum by the 
canonical momentum and finds the relation p = VS' — eA/c between the propagation 
and the wavefronts. In this case the wavefronts are no longer orthogonal to the classical 
trajectories. This is analogous to geometric optics in an anisotropic medium. 

The JWKB approximation is, in general, applicable when the change in the ampli¬ 
tude of the wave function is small at a scale of one wavelength, i.e. for a slowly varying 
amplitude of the wave function. This is usually not the case for reflections, or at classical 
turning points where fc —)■ 0. However, for most of these cases, methods were developed 
to calculate the additional phase shifts that are neglected when using the JWKB ap¬ 
proximation. Very good results are typically obtained by adding these additional phase 
shifts to the phase found by the JWKB approximation, even in cases where fc —)■ 0. 

1'2.2. Eikonal approximation. For most experiments in matter-wave optics the even 
simpler eikonal approximation of classical optics is sufficient. There, the phase of a wave 
function is calculated along the straight and unperturbed path between the starting point 
(source) and the observation point. 

1’3. Coherence. - Many of the phenomena in wave optics are concerned with the 
superposition of many waves. Therefore, one of the central questions is concerned with 
the coherence properties of this superposition. Naturally, this is also an important ques¬ 
tion in matter-wave optics. In general, one can define the coherence of matter waves 
analogously to the coherence in light optics, by using correlation functions. 

The first order correlation function for matter-waves with respect to coordinate rj and 
a displacement Sp is defined by 

(23) g^^HSr,) = {^\T^{dvM), 

where Tjj{Sr]) is the translation operator with respect to a displacement of drj. The width 
of this function with respect to the displacement rj is called the amount of coherence with 
respect to rj. 

1’3.I. Spatial coherence. In the case of spatial coherence, Tx(Jx) = exp(ip(5x//i) is 
the spatial translation operator and the correlation function takes the familiar form 

(24) (Sx) = {tp{x)\tp{x + 6x)) = J tp*{x)'tp{x + Sx) dx. 

(®) see for example chapter 6 in Messiah’s book on quantum mechanics [15] 
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In analogy to optics with light, in a beam of matter-waves one distinguishes between 
longitudinal and transverse coherence(®). 

Longitudinal coherence. An interesting and common example is the longitudinal coher¬ 
ence length Ic of a particle beam with a Gaussian distribution of fc-vectors, which prop¬ 
agates along the cc-direction. An example of such a beam is given in equation (9). Its 
first order longitudinal correlation function is given by 

6x^0-^ 

(25) g^^'^ (Sx) = {'tp{x)\tp{x + 6x)) = e 

and the longitudinal coherence length £c is related to the momentum distribution in the 
beam by 

(26) 4 = 1 = = (hd M 

(T k (T y 27r O'y 

where Oy is the rms velocity spread connected to the momentum distribution and (Ads) 
is the mean deBroglie wavelength. 

Because of the different dispersion relations, light and matter-waves exhibit differences 
in their correlation functions. For the linear dispersion relation of light propagating in 
vacuum, the coherence length can pictorially be associated with the size of a Fourier 
transform limited pulse with the same frequency width. For matter waves, the quadratic 
dispersion relation leads to a spreading of the wave packet, and the size of a wave packet 
can not be related to the coherence length, as illustrated in figure 1. This difference led 
to various discussions in the early matter-wave experiments with neutrons [16, 17, 18, 
19, 20]. 

It is interesting to note that the coherence length (transversal and longitudinal) has 
nothing to do with the size of the particles. For example for the matter-wave interferom¬ 
eter experiments with Na 2 molecules at MIT the coherence length was about a factor 4 
smaller than the size of the molecule (size of the Na 2 molecule ~ 400 pm, its de Broglie 
wavelength: Aae ~ 10 pm and its coherence length lyoh ~ 100 pm). Nevertheless, the 
same interference contrast as in the experiments with Na atoms was observed [21]. Sim¬ 
ilar conclusion at a more extreme scale can then be drawn from the later experiments on 
large molecules by the Vienna group (see the lectures of M. Arndt in these proceedings). 

Consequently a measurement of the longitudinal coherence length in a time-independent 
experiment using an interferometer [16, 17, 22] generally tells us nothing about the length 


(®) For massive particles the distinction between longitudinal and transversal coherence is not 
always as clear as for light. This can be easily seen if one notices that for a nonrelativistic 
particle longitudinal and transversal motion can be transformed into each other by a simple 
Galilean transformation. The distinction breaks down especially if the particles are brought to 
rest. Therefore, in the following discussion we will assume a particle beam with mean fc-vector 
k much larger than the momentum distribution (k ^ Uk). 
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Fig. 1. - Gaussian Wavepacket with fco = 10 and ak ~ 3 after a propagation time of t = 
10 (in natural units h = 1, m = 1). Att = 0 the wavepacket is a minimum uncertainty 
wavepacket with a width in x of — 0.3 {g). (a) Real part of the wave function Re(?/)(a;)). 
(b-e) Re('i/)*(a:)'i/)(a: + 5x)) for 4 different values of Sx. (b): Sx = 0. (c-e) for increasing Sx the 
oscillating contribution of Re(V'* {x)'4>{x + 5x)) leads to a vanishing contribution to the coherence 
function g^^^(Sx) displayed in (f). (g) Initial size of the wavepacket. 


or even the existence of a wave packet [18, 19, 20]. One can easily show that the longitudi¬ 
nal correlation functions for the following examples are identical: a minimum uncertainty 
Gaussian wave packet, the same wave packet after spreading for an arbitrary time, and 
even a superposition of plane waves with the same k-vector distribution but random 
phases tpr- The latter is the correct description for an thermal atomic source, as it is 
used in most atom optics experiments. Only in experiments when the beam is chopped at 
a timescale comparable to the inverse energy spread, can one hope to prepare an atomic 
beam in a not completely chaotic state. 


Transverse coherence. The transverse coherence of a matter-wave is obtained similar 
to the longitudinal coherence length, but by translation in the transverse direction 5z 
fftrans('^' 2 :) = {il}{z)\il}{z -\- 5z)). In analogy to light optics it is related to the transverse 
momentum distribution by the von Cittert-Zernike theorem [3]. For an atomic beam the 
transverse coherence of a beam can be defined by preparation (selection) in space (e.g. 
collimation by slits), transverse to the propagation direction. For waves emitted with an 
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angular spead a, the longitudinal coherence length is 


(27) 


^tcoh 


1 Ads 
2 ~ 


^dB Z 

2 w 


where z is the distance from the source aperture of width w. For a BEC, transverse 
coherence is related to fluctuations, and the resulting multimode structure of the quantum 
gas, similar to a laser that emits in multiple transverse modes. 

1’3.2. Coherence in momentum space. To address the problems arising in the inter¬ 
pretation of the experiments on the longitudinal coherence length for matter-waves a 
more useful concept is to study coherence in momentum space. In analogy to coherence 
in real space the correlation function in momentum space is given by 


(28) 


g^^\dk) 




J i;*{k)'tp{k + 6k)dk. 


This correlation function is, in principle, sensitive to the phase relations between different 
^-components of an atomic beam, and can therefore distinguish between the different 
interpretations of the longitudinal coherence length. 

This correlation function in momentum space can only be measured in a time-dependent 
experiment. One possibility the measure g^^'^{Sk) for matter-waves is the sideband inter¬ 
ferometer as described by B. Golub and S. Lamoreaux [23]. There, g^^'^{Sk) is measured 
by superposing two paths, where the energy of the propagating wave is shifted in one of 
the two(^). The result of g^^^Sk) is the given as the relative contrast of the time-averaged 
interference fringes in the interferometer. Similarly coherence in momentum space can be 
probed by a differentially tuned separated oscillatory field experiment as demonstrated 
at MIT [25] 

1’3.3. Higher-order coherence. The higher-order correlation functions for matter-waves 
can be defined similarly to light optics. The second-order correlation function is given 
by the joint detection probability of two particles at two locations: 

(29) g^'^\xi,ti;x2,t2) = {'S\b\xi,ti)b\x2,t2) b{xi,ti)b{x2,t2)\’i’) 

where b^ {x, t) and b{x, t) are the creation and annihilation operators for an atom at time 
t and location x, and is a multi-particle wave function. 

The second-order correlation function is trivially zero for one particle experiments. 
Since the form is similar to light optics one would expect bunching behavior for bosons 
and antibunching behavior for fermions from a chaotic source. High phase-space density 
per propagating mode is needed for realistic experiments to observe these higher-order 
coherences. 


(^) Experiments describing interference of neutron paths with different energies are described 
in ref. [24]. 
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Second- and higher-order correlation functions for identical particles emitted from 
highly-excited nuclei were investigated in refs. [26, 27] as a tool to measure the coherence 
properties of nuclear states. Second-order correlation functions for fermions were indeed 
observed in nuclear physics experiments by studying neutron correlations in the decay 
of highly-excited ^"^Ca nuclei [28]. Higher-order correlation functions in pions emitted in 
heavy-ion collisions were widely used to investigate the state of this form of matter [29, 
30]. In atom optics, the bunching of cold atoms in a slow atomic beam, analogous to the 
Hanbury-Brown-Twiss experiment in light optics [31, 32], was Hrst observed for Ne [33] 
and Li [34] and recently using ultracold "‘He [35]. It was measured for atom lasers [36] or 
to study the BEC phase transition [37, 38, 39] or higher-order correlation functions [40]. 
Fermi antibunching in an atomic beam of '^He was observed in ref. [41]. 

1’4. Index of refraction for matter waves. - We define the index of refraction for 
matter waves in the same way as for light: as the ratio of the free propagation fc-vector 
ko to the local fc-vector fc(r) 

(30) n(r) = 

fco 

There are two different phenomena that can give rise to a refractive index for matter- 
waves: 

• One can describe the action of a potential as being equivalent to a refractive index. 

• One finds a refractive index from the scattering of the propagating particles off a 
medium. This is equivalent to the refractive index in light optics. 

1’4.1. Index of refraction caused by a classical potential. If we compare the local fc- 
vector in equation (5) to our above definition of a refractive index (equation (30)), one 
sees that in matter-wave optics the action of a scalar potential V (r) can be described 
as a position dependent refractive index n(r). In most cases the potential Vir) will be 
much smaller than the kinetic energy Ekin (b^(r) ^ Ekin) of the atom, and one finds 


(31) 


n(r) 



^(r) 


Therefore, the refractive index will be larger than unity (n(r) > 1) in regions with 
an attractive potential (H(r) < 0) and smaller than unity (n(r) < 1) in regions with a 
repulsive potential. It is interesting to note that the refractive index caused by a classical 
potential has a strong dispersion, as n changes with (1/fco)- 

The above relation in equation (31) can be extended to a vector potential A. In this 
case, the canonical momentum Pc (wave vector kc) and the kinetic momentum pk (wave 
vector kfc) are not necessary parallel. In a wave description, the wavefront (orthogonal 
to kc) and the propagation direction (parallel to k^) are not orthogonal to each other. 
This is similar to propagation of light in an anisotropic medium. The refractive index is 
then direction dependent. 
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1'4.2. Index of refraction from scattering. A second phenomena that gives rise to a 
refractive index for matter-waves is the interaction of the matter-wave with a medium. 
This is analogous to the refractive index for light, which results from the coherent forward 
scattering of the light in the medium. Similarly, the scattering processes of massive 
particles inside a medium result in a phase-shift of the forward scattered wave and define 
a refractive index for matter-waves. Here we give a schematic introduction. For a full 
treatment, see one of the standard books on scattering theory. 

From the perspective of wave optics the evolution of the wave function i/i while prop¬ 
agating through a medium is given by 

(32) i^(x) = ■i/l(0) gi If Wa;Re(/(fcc,0)) g-f)) Aa: Im(/(fc„ .0)) ^ 

Here /ciab is the wave vector in the laboratory frame, kc the wave vector in the center-of- 
mass frame of the collision, N is the areal density of scatterers in the medium and f{kc,0) 
is the center-of-mass, forward scattering amplitude. The amplitude of propagating wave 
function ip is attenuated in proportion to the imaginary part of the forward scattering 
amplitude, which is related to the total scattering cross section by the optical theorem 

Air 

(33) (Ttot =-^lm(/(fcc,0)). 

Kc 

In addition, there is a phase shift (p proportional to the real part of the forward scattering 
amplitude 

O'jr 

(34) (p{x) = — NxRe{f{kc,0))). 

Kc 

In analogy to light optics one can define the complex index of refraction 

Ott 

(35) n = l + - - —Nf{kc,0). 

Mab 

The refractive index of matter for de Broglie waves has been extensively studied in 
neutron optics [42, 43], especially using neutron interferometers. It has also been widely 
used in electron holography [44]. In neutron optics, scattering is dominantly s-wave and 
measuring the refractive index gives information about the s-wave scattering length a 
defined as 

(36) a = - lim /o, 

k—^0 

where /o is the s-wave scattering amplitude. 

In atom optics with thermal beams, usually many partial waves, typically a few 
hundred, contribute to scattering of thermal atoms. The number of contributing partial 
waves I can be estimated by I ~ where Xj. is the range of the inter atomic potential 
and kc is the center-of-mass wave vector of the collision. The refractive index will depend 
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on the forward scattering amplitude and therefore on details of the scattering process. 
Measuring n will lead to new information about atomic and molecular scattering [45], 
especially the real part of the scattering amplitude, not directly accessible in standard 
scattering experiments. 

For ultracold atoms and an ultracold media the scattering is predominately s-wave 
and can be described by a scattering length a very similar to neutron optics. This 
regime can be reached for scattering inside a sample of ultracold atoms like a BEC 
or by scattering between two samples of ultracold atoms. We would like to note that 
for scattering processes between identical atoms at ultra low energies, quantum statistic 
becomes important. For example scattering between identical Fermions vanishes because 
symmetry leads to suppression of s-wave scattering. On the other hand the dominance 
of s-wave scattering at low energies is only valid if the interaction potential decays faster 
than 1/r^. For scattering of two dipoles, higher partial waves contribute even in the limit 
of zero collision energy. 

We now look closer at the low-energy limit where s-wave scattering is predominant. 
The scattering can in first approximation be described by only one parameter, the scat¬ 
tering length a. Here we can derive simple relations for the dispersion of the refractive 
index n starting from a low energy expansion of the s-wave scattering amplitude 

(37) f{k,0) = - 1) ~ -a(l - ika) 

2ik 

where Jq — —ka is the s-wave scattering phase-shift(®). The refractive index is then 
given by 

27r 

(38) n(fciab) = 1 + ^- —Na{l-ika) 

Mab 

For fc —)■ 0 n becomes predominantly real and diverges with 1/fc^. 

Consequently we can now reverse the above argument defining a refractive index for a 
classical potential and define an effective optical potential C/opt for a particle in a medium 
with scattering length a 

(39) Copt(r) = 

nic 

where rric is the reduced mass. This potential C/opt is one of the basics ingredients for 
many neutron optics experiments and neutron optics devices. 

An important phenomenon in matter-wave optics, is that matter interacts with itself, 
matter-wave optics is inherently non-linear, and the non-linearities can define the dom¬ 
inant energy scale. The local refractive index, and therefore the propagation of matter 


(®) To be more precise for larger k one can use the effective range approximation fccot((5o) = 
— 1/a -I- (refi is the effective range of the potential 
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waves depends on the local density of the propagating particles. A simple description 
can be found in the limit when the propagating beam can be viewed as weakly inter¬ 
acting, that is if the mean particle spacing is much larger than a {a where p is 

the density). The self-interaction can then be described by the optical potential (equa¬ 
tion (39)). In its simplest form, it leads to an additional term in the Schrodinger equation 
(equation (3)) which is then nonlinear and called the Gross-Pitaevskii equation [46] 


(40) 



V{r,t) 




a 

m 


= ih 


dtp{r, t) 
dt 


This self-interaction leads to a new type of nonlinear optics where even the vacuum is 
nonlinear. This has to be contrasted with the fact that for light, nonlinearities are very 
small and come into play only in special media. 
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2. Optics and interferometry using gratings 

In this section we will give an overview of optics and interferometry with beams of 
atoms or molecules using gratings. We will only discuss the main aspects and phenom¬ 
ena, and refer the reader for details about experiments to the review article on atom 
interferometry by Cronin et al. [1]. 

2’1. Diffraction. — Diffraction of matter waves from phase and amplitude modulating 
objects is a hallmark example of wave propagation and interference. It arises from the 
coherent superposition and interference of the propagating matter wave which is modified 
in amplitude and phase by the diffracting structure. It is described by the solution of 
the Schrodinger equation (equation (3)) with the appropriate boundary conditions. An 
elegant approach to solve this problem is to express the amplitude and phase of the 
matter wave at a position in space and time as the sum over all possible paths between 
the source and the observation point (see section 1'2). Beamsplitters for atom beam 
interferometers are often based on diffraction. Comparing the diffraction of matter waves 
and light, we expect differences arising from the different dispersion relations. These 
manifest themselves in time dependent diffraction problems, and will give rise to a new 
phenomenon: diffraction in time. 

2T.1. Diffraction in space. First we discuss diffraction in space, transverse to the 
propagation of the beam. A diffraction grating is a diffracting region that is periodic. 
Spatial modulation of the wave by the grating generates multiple momentum components 
for the scattered waves which interfere. The fundamental relationship between the mean 
momentum transferred to waves in the component and the grating period, d, is 

(41) 6pn = n— = nhG 6n ~ 

where G = 2TT/d is the reciprocal lattice vector of the grating, h is Planck’s constant, 
and AdB = h/pbeam is the de Broglie wavelength of the incoming beam. In the far field 
diffraction is observed with respect to the diffraction angle 6*„. To resolve the different 
diffraction orders in the far field the transverse momentum distribution of the incoming 
beam must be smaller than the transverse momentum given by the diffraction grating 
6p = hG. This is equivalent to the requirement that the transverse coherence length 
must be larger than a few grating periods. This is usually accomplished by collimating 
the incident beam(®). 

The first examples of atom interference were diffraction experiments. Just three 
years after the electron diffraction experiment by Davisson and Germer [4] Estermann 


(®) The transverse coherence length is Itcoh ~ AdB/i?coii, where Arb is the de Broglie wavelength 
and i?coii is the (local) collimation angle of the beam (the angle subtended by a collimating slit). 
Since for thermal atomic beams Arb ~ 10 pm a collimation of dcoii < lOprad is required for a 
1 pm coherent illumination. 
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Fig. 2. - (Left) Historic data showing diffraction of He atoms from a LiF crystal surface [8]. 
(Right) Diffraction of He atoms transmitted through a nanostructure grating. The average 
velocity, the velocity spread of the beam, the uniformity of the material grating, and the strength 
of atom-sruface van der Waals forces can all be determined from these data [47]. The insert 
shows a 100 nm period grating for atom waves. The dark regions are slits, and light regions 
are free-standing silicon nitride bars. Figure courtesy of T.A. Savas and H.I. Smith at the MIT 
Nanostructure laboratory [48, 49, 50]. Figure curtesy of J.P. Toennies, W. Schoellkopfand and 
O. Kornilov. 


and Stern observed diffraction of He beam off a LiF crystal surface [8]. 

Classical wave optics recognizes two limiting cases, near- and far-field. In the far- 
field the curvature of the atom wave fronts is negligible and Fraunhofer diffraction is a 
good description. The diffraction pattern is then given by the Fourier transform of the 
transmission function, including the imprinted phase shifts. In the near-field limit the 
curvature of the wave fronts must be considered and the intensity pattern of the beam 
is characterized by Fresnel diffraction. Edge diffraction and the Talbot self-imaging of 
periodic structures are typical examples. 

2T.2. Diffraction from nano-fabricated structures. With the advent of modern nan¬ 
otechnology it became possible to fabricate elaborate arrays of holes and slits with feature 
sizes well below 50 nm in a thin membrane that allow atoms to pass through. Diffraction 
from a fabricated grating was first observed for neutrons by H. Kurz and H. Rauch in 
1969 [51] and for atoms by the Pritchard group at MIT [52]. The latter experiment 
used a transmission grating with 200 nm wide slits. Similar mechanical structures have 
been used for single slits, double slits, diffraction gratings, zone plates, hologram masks, 
mirrors, and phase shifting elements for atoms and molecules. 

The benefits of using mechanical structures for atom optics include the possibility to 
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Fig. 3. - (Left) A zone plate for focusing atom beams. The plate has free-standing annular rings 
and radial support struts. (Right) Data showing focused and defocused atom beam components. 
Figure from ref. [53]. (Left) Image reconstruction by atom holography. (Middle Top) The 
hologram is designed by computer and realized with a SiN membrane with square holes. (Right) 
data from a switchable hologram. [54, 55, 56]. 


create feature sizes smaller than light wavelengths, arbitrary patterns, rugged designs, 
and the ability to diffract any atom or molecule. The primary disadvantage is that atoms 
or molecules can stick to (or bounce back from) surfaces, so that most structures serve 
as absorptive atom optics with a corresponding loss of transmission. When calculating 
the diffraction patterns, one has to consider that, first, nanofabrication is never perfect, 
and that the slits and holes can have variations in their size. Second, the van der Waals 
interaction between atoms and molecules and the material of the gratings can lead to 
effectively much smaller slits and holes in the diffracting structures. Moreover, the wave- 
front emerging from the hole can have additional phase shifts from the van der Waals 
interaction with the surface. Such effects can be particularly significant for molecules 
with a large electric polarizability and very small diffraction structures. For a detailed 
discussion of this topic we refer the reader to the lecture of M. Arndt. 

Nano-fabricated gratings have been used to diffract atoms and molecules such as ^He, 
^He 2 , ^Hea and larger ^He clusters (Fig. 2), Na 2 , Ceo, C 60 F 48 , C 44 H 30 N 4 and many 
more. [21, 57, 58, 59, 60, 61, 62, 63]. Fresnel zone plates have been employed to focus 
atoms [64, 53] and spot sizes below 2 /im have been achieved. Atom holography with 
nanostructures can make the far-field atom flux resemble arbitrary patterns. Adding 
electrodes to a structure allows electric and magnetic fields that cause adjustable phase 
shifts for the transmitted atom waves. With this technique, a two-state atom holographic 
structure was demonstrated [65, 66 , 56] that produced images of the letters </> or tt as 
shown in Fig. 3. The different holographic diffraction patterns are generated depending 
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Fig. 4. - Implementing complex potentials with light: (a) Top: Schematic of an open transi¬ 
tion from a metastable state. If the ground state is not detected then emission of a photon 
corresponds to ’absorption’ of the atom, and one can realize also ’absorptive’ potentials. Bot¬ 
tom-. real and imaginary part of the optical potential close to a resonance as depicted above, 
(b) Level structure of Ar atom with relevant transitions as an example to implement complex 
optical potentials. 


on the voltages applied to each nanoscale aperture. 

2’1.3. Light gratings from standig waves. In an open two-level system the interaction 
between an atom and the light field (with detuning A = wiaser — Watom) can be described 
by an effective optical potential of the form [67] (figure 4): 


(42) 


U{x) a I[x) 


1 

2A-i-ir 


where F is the atomic decay rate and I{x) is the light intensity. The imaginary part of 
the potential results from the spontaneous scattering processes, the real part from the ac 
Stark shift. If the spontaneous decay follows a path to a state which is not detected, the 
imaginary part of the potential in equation (42) is equivalent to absorption. On-resonant 
light can therefore be used to create absorptive structures. Light with large detuning 
produces a real potential and therefore acts as pure phase object. Near-resonant light 
can have both roles. 

The spatial shape of the potential is given by the local light intensity pattern, /(cc), 
which can be shaped with all the tricks of near and far field optics for light, including 
holography. The simplest object is a periodic potential created by two beams of light 
whose interference forms a standing wave with reciprocal lattice vector G = fci — /c 2 . 
Such a periodic light field is often called light crystal or more recently an optical lattice 
because of the close relation of the periodic potentials in solid state crystals, and thus 
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Fig. 5. - (a) Kapitza Dirac (KD) diffraction (b) Bragg Diffraction. The top row shows the 
essential difference: thick vs. thin gratings. The bottom row shows data obtained by the 
Pritchard group for KD and Bragg diffraction [68, 69]. 


motivates the use of Bloch states to understand atom diffraction. 

Since light gratings can fill space, they can function as either thin or thick optical ele¬ 
ments. In the case of a grating, the relevant scale is the Talbot length Txaibot = d^/Ade- 
For a spatial extent D along the beam propagation of D Txaibot it is considered thin, 
otherwise thick (figure 5). For a thin optical element the extent of the grating along the 
propagation direction has no influence on the final diffraction (interference). This limit is 
also called the Raman-Nath limit. For a thick optical element the full propagation of the 
wave throughout the diffracting structure must be considered. A thick grating acts as a 
crystal, and the characteristics observed are Bragg scattering or channeling, depending 
on the height of the potentials. 

The second distinction, relevant for thick gratings, has to do with the strength of the 
potential. One must determine if the potential is only a perturbation, or if the potential 
modulations are larger then the typical transverse energy scale of the atomic beam or 
the characteristic energy scale of the grating, 

(43) Eg = h^G^/(2m) = Ahujrec, 

associated with one grating momentum unit HG (huj-ec is an atom’s ‘recoil energy’ due 
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Fig. 6. - (color online) Diffraction from a measurement-induced grating, (a) Schematic of two 
on-resonant standing waves of light. The first causes atom diffraction. The second translates 
to analyze near-held atomic hux. (b) Periodic structure in the transmitted atomic beam, (c) 
Far-held atom diffraction from a measurement induced grating. Figure from ref. [70]. 


to absorbing or emitting a photon). For weak potentials, U Eq, one observes Bragg 
scattering. The dispersion relation looks like that of a free particle, only with avoided 
crossings at the edges of the cell boundaries. Strong potentials, with U Eq, cause 
channelling. The dispersion relations are nearly flat, and atoms are tightly bound to the 
wells. 


Diffraction with on-resonant light. If the spontaneous decay of the excited state proceeds 
mainly to an internal state which is not detected, then tuning the light frequency of 
a standing light wave to resonance with an atomic transition (A = 0) can make an 
‘absorptive’ grating with light. (If the excited state decays back to the ground state, 
this process produces decoherence.) For a thin standing wave the atomic transmission is 
given by 


(44) 


T{x) =exp 


-^[l-bcos(Gx)] , 


where k is the absorption depth for atoms passing through the antinodes. For sufficiently 
large absorption only atoms passing near the intensity nodes survive in their original state 
and the atom density evolves into a comb of narrow peaks. Since the ‘absorption’ in¬ 
volves spontaneous emission, such light structures have been called measurement induced 
gratings. As with all thin gratings, the diffraction pattern is then given by the scaled 
Fourier transform of the transmission function. 

On-resonant standing waves have been used as gratings for a series of near-held (atom 
lithography, Talbot effect) and far-held (diffraction, interferometry) experiments [70, 71, 
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72, 73, 74]. An example is shown in figure 6. These experiments demonstrate that 
transmission of atoms through the nodes of the ‘absorptive’ on resonant light masks is a 
coherent process. 

Light crystals. If the standing wave is thick, one must consider the full propagation of the 
matter wave inside the periodic potential. The physics is characterized by multi-wave 
(beam) interference. For two limiting cases one can get to simple models. For weak 
potentials: Bragg scattering; and for strong potentials: coherent channeling. 

When an atomic matter wave impinges on a thick but weak light crystal, diffraction 
occurs only at specific angles, the Bragg angles 9b defined by the Bragg condition 

(45) A^AdB = Ap,,sin(6»B). 

Bragg scattering, as shown in figure 5 (right column) transfers atoms with momentum 
—Px into a state with a single new momentum, px = —px + hG. Momentum states in 
this case are defined in the frame of the standing wave in direct analogy to neutron 
scattering from perfect crystals. Bragg scattering of atoms on a standing light wave 
was first observed at MIT [69]. Higher order Bragg pulses transfer multiples of NhG of 
momentum, and this has been demonstrated up to 50*^ order (transter of IQOhk and 
beyond [75, 76, 77, 78, 79] (see also cectures by G. Tio, M. Kasevich and H. Mueller) 

Bragg scattering can be described as a multi-beam interference as treated in the 
dynamical diffraction theory developed for neutron scattering. Inside the crystal one 
has two waves, the incident ‘forward’ wave (kp) and the diffracted ‘Bragg’ wave (fcs). 
These form a standing atomic wave with a periodicity that is the same as the standing 
light wave. This is enforced by the diffraction condition (fcs — kp = G). At any location 
inside the lattice, the exact location of atomic probability density depends on kp, ks 
and the phase difference between these two waves. 

For incidence exactly on the Bragg condition the nodal planes of the two wave fields 
are parallel to the lattice planes. The eigenstates of the atomic wave field in the light 
crystal are the two Bloch states, one exhibiting maximal (ipmax) the other minimal (ppmin) 
interaction: 


'^max 

(46) ’4^min 

For tpmax the anti-nodes of the atomic wave field coincide with the planes of maximal light 
intensity, for ipmin the anti-nodes of atomic wave fields are at the nodes of the standing 
light wave. These states are very closely related to the coupled and non-coupled states 
in velocity selective coherent population trapping (VSCPT). 

The total wave function is the superposition of ipmax and ipinin which satisfies the 
initial boundary condition. The incoming wave is projected onto the two Bloch states 


I r _ j 

— g''2 _|_g ‘'2 
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mirror angle [grad] maximum of light field 

Fig. 7. - (color online) Bragg diffraction of atoms from resonant standing waves of light, (a) 
Anomalous transmission of atoms, (b) A resonant standing wave inside an off-resonant light 
crystal serves to measure the atom wave fields inside the light crystal. Figure from ref. [67]. 


which propagate through the crystal accumulating a relative phase shift. At the exit, 
the final populations in the two beams is determined by interference between ij^max and 
ipmin and depends on their relative phase. 

Bragg scattering can also be observed with absorptive, on-resonant light structures [67] 
and combinations of both on- and off-resonant light fields [80]. One remarkable phe¬ 
nomenon is that one observes the total number of atoms transmitted through a weak 
on-resonant standing light wave increases if the incident angles fulfill the Bragg condition, 
as shown in figure 7. This observation is similar to what Bormann discovered for X-rays 
and called anomalous transmission [81]. It can easily be understood in the framework 
of the two beam approximation outlined above. The rate of de-population of the atomic 
state is proportional to coupling between the atom wave field and the standing light field. 
The minimally coupled state ijjmin will be attenuate much less and propagate much fur¬ 
ther into the crystal than 'll)max- For a sufficiently thick light crystal the propagating 
wave field will be nearly pure V’mm at the exit, and as a consequence one observes two 
output beams of equal intensity. Inserting an absorptive mask [70] inside the light crystal 
one can directly observe the standing matter wave pattern inside the crystal [67, 82] and 
verify the relative positions between the light field and ipmin- 

More complex potentials for atoms can be made by superposing different standing 
light waves [80]. For example, a superposition of a on- and a off-resonant standing wave 
with a phase shift Aip = ±7r/2 results in a combined potential of U{x) = which, 

in contrast to a standing wave, has only one momentum component and therefore only 
diffract in one direction [80]. 

2T.4. Diffraction in time. If the optical elements are explicitly time dependent an 
interesting phenomenon arises for matter waves, which is generally not observed in light 
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Fig. 8. - The ENS diffraction in time experiments described in the text. Figure from ref. [88] 
and reviewed in ref. [89]. 


optics: diffraction in time. The physics behind this difference is the different dispersion 
relations. The time dependent optical element creates different energies which propagates 
differently according the quadratic dispersion relation for matter waves, which leads to 
superposition between matter waves emitted at different times at the detector. 

Diffraction in time and the differences between light and matter waves can best seen 
by looking at the shutter problem as discussed by Moshinsky [83]. We start with a 
shutter illuminated with monochromatic wave from one side. The shutter opens at time 
t = 0 and we ask the question how does the time dependence of the transmitted radiation 
behave at a distance 2 ; behind the shutter. For light we expect a sharp increase in the light 
intensity at time t = zfc. For matter waves the detected intensity will increase at time 
t = z/v, where v is the velocity corresponding to the incident monochromatic wave. But 
the increase will be not instantaneous, but will show a typical Fresnell diffraction pattern 
in time, similar to the diffraction pattern obtained by diffracting from an edge [83, 84]. 
Similar arguments also hold for diffraction from a single slit in time, a double slit in time, 
or any time structure imprinted on a particle beam. 

The experimental difficulty in seeing diffraction in time is that the time scale for 
switching has to be faster than the inverse frequency (energy) width of the incident 
matter wave. This condition is the time equivalent to coherent illumination in spatial 
diffraction. The first (explicit) experiments demonstrating diffraction in time used ultra¬ 
cold neutrons reflecting from vibrating mirrors [85, 86, 87]. Side bands of the momentum 
components were observed. 

The group of J. Dalibard at the ENS in Paris used ultra cold Cs atoms (T ~ 3.6 /iK) 
released from an optical molasses reflecting from an evanescent wave atom mirror [88, 
90, 91]. By pulsing the evanescent light field one can switch the mirror on and off, 
creating time-dependent apertures that are diffractive structures. To obtain the necessary 
temporal coherence a very narrow energy window was first selected by two (0.4 ms) 
temporal slits separated by 26 ms. If the second slit is very narrow (< 10 /is) one observes 
single slit diffraction in time, if the mirror is pulsed on twice within the coherence time 
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of the atomic ensemble one observes double slit interference in time, and many pulses 
lead to a time-dependent flux analogous to a grating diffraction pattern as shown in 
figure 8. From the measurement of the arrival times of the atoms at the final screen by 
fluorescence with a light sheet the energy distribution can be reconstructed. 

Because the interaction time between the atoms and the mirror potential (< 1 /is) 
was always much smaller then the modulation time scale (> 10 ^s), these experiments 
are in the ‘thin grating’ (Raman-Nath) regime for diffraction in time. 


Modulated light crystals. The time equivalent of spatial Bragg scattering can be reached if 
the interaction time between the atoms and the potential is long enough to accommodate 
many cycles of modulation. When a light crystal is modulated much faster than the 
transit time, momentum is transferred in reciprocal lattice vector units and energy in 
sidebands at the modulation frequency. This leads to new resonance conditions and 
’Bragg diffraction’ at two new incident angles [92, 93]. Consequently Bragg scattering in 
time can be understood as a transition between two energy and momentum states. The 
intensity modulation frequency of the standing light wave compensates the detuning of 
the Bragg angle and frequency of the de Broglie wave diffracted at the new Bragg angles 
is shifted by zLHuJmod [92, 93]. Thus, an amplitude modulated light crystal realizes a 
coherent frequency shifter for a continuous atomic beam. It acts on matter waves in an 
analogous way as an acousto-optic modulator acts on photons, shifting the frequency 
(kinetic energy) and requiring an accompanying momentum (direction) change. In a 
complementary point of view the new Bragg angles can be understood from looking at 
the light crystal itself. The modulation creates side bands ±Wmod on the laser light, and 
creates moving crystals which come from the interference between the carrier and the 
side bands. Bragg diffraction off the moving crystals occurs where the Bragg condition 
is fulfilled in the frame co-moving with the crystal,resulting in diffraction of the incident 
beam to new incident angles. 

The coherent frequency shift of the Bragg diffracted atoms can be measured by in¬ 
terferometric superposition with the transmitted beam. Directly behind the light crystal 
the two outgoing beams form an atomic interference pattern which can be probed by a 
thin absorptive light grating [70]. Since the energy of the diffracted atoms is shifted by 
fiwmod 7 the atomic interference pattern continuously moves. This results in a temporally 
oscillating atomic transmission through the absorption grating (see figure 9). 

Starting from this basic principle of frequency shifting by diffraction from a time de¬ 
pendent light crystal many other time-dependent interference phenomena were studied 
for matter waves [94, 93] developing a diffractive matter wave optics in time. For ex¬ 
ample using light from two different lasers one can create two coinciding light crystals. 
Combining real and imaginary potentials can produce a driving potential of the form 
U{t) ~ -which contains only positive (negative) frequency components respec¬ 

tively. Such a modulation can only drive transitions up in energy (or down in energy). 

Figure 10 summarizes thick and thin gratings in space and also in time with Ewald 
constructions to denote energy and momentum of the diffracted and incident atom waves. 
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c) Matter wave beating 
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Fig. 9. - (color online) Frequency shifter for matter waves, (a) A time-modulated light crys¬ 
tal causes diffraction in time and space, (b) Rocking curves show how the Bragg angle for 
frequency-shifted matter waves is controlled by the grating modulation frequency, (c) Matter 
wave interference in time. Figure from ref. [92]. 


The diffraction from (modulated) standing waves of light can also be summarized with 
the Bloch band spectroscopy picture [93, 95]. 

2’2. Interferometers. — Interferometers are, very generally speaking, devices that uti¬ 
lize the superposition principle for waves and allow a measurement through the resulting 
interference pattern. A generic interferometer splits an incoming beam ['0) in (at least 
two) different components \ipi), |' 02 )v which evolve along different paths in configura¬ 
tion space and then recombined to interfere. Interferometer exhibit a closed path and 
can be viewed, in a topological sense, as a ring. 

At the output port of an interferometer the superposition of two interfering waves 
V'out = V'l + V '2 leads to interference fringes in the detected intensity: 

I = = IV’l + V’2p 

(47) = Af -f A2 2A1A2 cos{(j)) 

= (/)(l + C7cos(</.)) 


where Ai = [V'lP, A 2 = |' 02 p are the amplitudes of the interfering beams and (j) = ipi—ip 2 
is their phase difference. The observed interference is then fully characterized by its phase 
Aip, and by two of the following: its amplitude 2 A 1 A 2 , its average intensity (/) = Af+A^, 
or its contrast C given by: 


(48) 


C = 


2A1A2 

Af + Al 
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Fig. 10. - Momentum diagrams for cases: (A) a thick grating, (B) a thin grating, (C) a thick 
pulsed grating (D) a thick harmonically modulated grating. Figure from ref. [93]. 


If one of the interfering beams is much stronger then the other, for example 
then the contrast of the interference pattern scales like 


(49) 



Consequently one can observe 20% (2%) contrast for an intensity ratio of 100:1 (10^:1) 
in the interfering beams. 

While the overall phase of the wave function is not observable, the power of inter¬ 
ferometers lies in the possibility of measuring the phase difference Atp between waves 
propagating along two different paths. 

(50) = (^1 - (^2 = ^{Si - S 2 ) 

h 


where (pi (Si) and <^2 {S 2 ) are the phases (classical action) along path 1 (path 2) of 
the interferometer. The other parameters of the interference pattern, the amplitude 
2 A 1 A 2 , the contrast C, and the mean intensity (I) also give information about the 
paths. Especially the contrast of the interference pattern tells us about the coherence in 
the interferometer. 

From a practical point of view interferometers can be divided in two categories: 

• In internal state interferometers, the beam splitter produces a superposition of in¬ 
ternal states which can be linked to external momentum. Examples are polarisation 
interferometry with light and Ramsey spectroscopy for internal states of massive 
particles. 
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• In de Broglie wave interferometers the beam splitter does not change the internal 
state but directly creates a superposition of external center of mass states and thus 
distinctly different paths in real space. The spatially distinct interfering paths can 
be created by wavefront division like in Young’s double slit or by amplitude division 
as realizes by a beam splitter in optics or by diffraction. In matter wave optics 
amplitude division Mach-Zehnder interferometer can be built with three gratings. 

When designing and building interferometers for beams of atoms and molecules, one 
must consider their specifics. (1) Beams of atoms and molecules have a wide energy 
distribution and consequently the coherence lengths for matter waves is in general very 
short (~I 00 pm for thermal atomic beams, and seldom larger then 10 fj,m even for 
atom lasers or BEC). This requires that the period and the position of the interference 
fringes must be independent of the de Broglie wavelength of the incident atoms. In 
optical parlance this is a property of white light interferometers like the three grating 
Mach-Zehnder configuration. (2) Atoms interact strongly with each other. Therefore 
the optics with matter waves is non-linear, especially in the cases where the atoms have 
significant density as in a BEC or atom laser. (3) Atoms can be trapped which allows 
a different class of interferometers for confined particles, which will be disused in a later 
section. 

2’2.1. Three-grating Mach-Zehnder interferometer. The challenge of building a white 
light interferometer for matter waves is most frequently met by the 3-grating Mach 
Zehnder (MZ) layout. In symmetric setup the fringe spacing is independent of wavelength 
and the fringe position is independent of the incoming direction (^*^). This design was 
used for the first electron interferometer [96], for the first neutron interferometer by H. 
Rauch [5], and for the first atom interferometer that spatially separated the atoms [ 6 ]. In 
the 3-grating Mach Zehnder interferometer the role of splitter, recombiner and mirror is 
taken up by diffraction gratings. At the position of the third grating (G3) an interference 
pattern is formed with the phase given by 

(51) (j) = G{xi - 2X2 + Xs) + Acfint 

where xi, X 2 , and X 3 are the relative positions of gratings 1, 2 and 3 with respect to an 
inertial frame of reference [97]. 

It is interesting to note that many diffraction-based interferometers produce fringes 
even when illuminated with a source whose transverse coherence length is much less than 
the (large) physical width of the beam. The transverse coherence can even be smaller 
then the grating period. Under the latter condition, the different diffraction orders will 


Diffraction separates the split states by the lattice momentum, then reverses this momentum 
difference prior to recombination. Faster atoms will diffract to smaller angles resulting in less 
transverse separation downstream, but will produce the same size fringes upon recombining with 
their smaller angle due to their shorter deBroglie wavelength. For three evenly spaced gratings, 
the fringe phase is independent of incident wavelength. 
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Fig. 11. - (a) and (b) Atom Interferometer setup used in Keith et al. [6]. (c) Interference fringe 
data and best fit with (!) = 157,000 counts per second and C = 0.42. A total of 5 seconds of 
data are shown and the uncertainty in phase is at/, = 2.7 x 10“® radians, (d) Average intensity 
(/) and contrast C as a function of detector position [(under different conditions than (c)]. 


not be separated and the arms of the interferometer it will specially overlap. Nevertheless, 
high contrast fringes will still be formed. 

The three grating interferometer produces a “position echo” as discussed by CDK [98]. 
Starting at one grating opening, one arm evolves laterally with hG more momentum 
for some time, the momenta are reversed, and the other arm evolves with the same 
momentum excess for the same time, coming back together with the other arm at the 
third grating. If the gratings are registered, its trapezoidal pattern starts at a slit on the 
first grating, is centered on either a middle grating slit or groove, and recombines in a 
slit at the third grating. Not surprisingly, spin-echo and time-domain echo techniques 
(discussed below) also offer possibilities for building an interferometer that works even 
with a distribution of incident transverse atomic momenta. 

Interferometer with nano fabricated gratings. The first 3-grating Mach-Zehnder interfer¬ 
ometer for atoms was built by Keith et al. [6] using three 0.4-/im period nano fabricated 
diffraction gratings. Starting from a supersonic Na source with a brightness of i? « 10^® 
s“^cm“^sr“^ the average count rate (/), in the interference pattern was 300 atoms per 
second. Since then, gratings of 100 nm period have been used to generate fringes with 
300000 atoms per seconds. 

Following the design shown in Fig. 11) there are two the MZ Interferometers formed 
starting with a common incident beam. One by two paths created by O**', -1-1’**^ and 
-1-1®*, —1®* order diffraction at grating Gl, G2 respectively. The second one formed 
symmetrically by two paths created by 0**“, —1®* and —1®*, -1-1®* order diffraction at 
grating Gl, G2 respectively. In each interferometer loop the difference in momentum is 
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one unit of HGx and at the position of the third grating G3 an interference pattern forms 
in space as a standing matter wave wave with a period oi d = G/{2 t:) and a phase that 
depends on the location xi and X 2 of the two gratings G1 and G2 as well as the interaction 
phase A<f>int- These fringes can be read out in different ways. The simplest is to use 
the third grating as a mask to transmit (or block) the spatially structured matter wave 
intensity. By translating G3 along x one obtains a moire filtered interference pattern 
which is also sinusoidal and has a mean intensity and contrast 


(52) 


^ _ sin(Gu;3/2) ~ 
d {GW 3 / 2 ) ^ 


where / and G refer to the intensity and contrast just prior to the mask. 

There are in fact many more interferometers formed by the diffraction from absorption 
gratings. For example, the and 2"'^ orders can recombine in a skew diamond to pro¬ 
duce another interferometer with the white fringe property. The mirror images of these 
interferometers makes contrast peaks on either side of the original beam axis (figure 11). 
In the symmetric MZ interferometer all those interferometers have fringes with the same 
phase, and consequently one can therefore build interferometers with wide uncollimated 
beams which have high count rate, but lower contrast (the contrast is reduced because 
additional beam components which do not contribute to the interference patterns such 
as the zeroth order transmission through each grating will also be detected). 

For well-collimated incoming beams, the interfering paths can be separated at the 2"'^ 
grating. For example in the interferometer built at MIT the beams at the 2”^^ grating 
have widths of 30 /im and can be separated by 100 /rm (using 100-nm period gratings and 
1000 m/s sodium atoms (Ade = 16 pm). Details of this apparatus, including the auxiliary 
laser interferometer used for alignment and the requirements for vibration isolation, are 
given in ref. [97]. 


Interferometers with light gratings. One can also build MZ interferometers with near¬ 
resonant standing light waves which make species-specific phase gratings (Fig. 12). The 
third grating can function to recombine atom waves so their relative phase dictates 
the probability to Hnd atoms in one output port (beam) or another. Alternatively, 
fringes in position space can be detected with fluorescence from a resonant standing 
wave. Another detection scheme uses backward Bragg scattering of laser light from the 
density fringes. Detecting the direction of exiting beams requires that the incident beams 
must be collimated well enough to resolve diffraction, and may well ensure that the beams 
are spatially separated in the interferometer. Because they transmit all the atoms, light 
gratings are more efficient than material gratings. 

Rasel et al. [74] used light gratings in the Kapitza-Dirac regime with a 5/xm wide 
collimated beam. Many different interferometers are formed, due to symmetric KD 
diffraction into the many orders. Two slits after the interferometer served to select both 
the specific interferometer, and the momentum of the outgoing beam (ports 1 and 2 in 
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Fig. 12. - Atom interferometers based on three standing waves of light, (a) Atom beam and 
three Kapitza-Dirac gratings, (b) Atominterference patterns for both output ports demonstrate 
complementary intensity variations. This is a consequence of atom number conservation. Figures 
a and b reproduced from [74]. (c) Interferometer based on three Bragg gratings. Dashed 

line shows the path of auxiliary optical interferometer used for stabilization, (d) Intensity 
fluctuations in beam A vs. position of the Bragg gratings. For second order Bragg diffraction, 
fringes of half the period are formed. Figures c and d reproduced from [99]. (e) Schematic of 
the 7r/2 — tt — 7r/2 Bragg interferometer for atoms in a BEC falling from a trap. (f)Absorption 
images and density prohles demonstrating different outputs of the interferometer. Figures e and 
f reproduced from [100]. 


figure 12a). Fringes show complementary intensity variations, as expected from particle 
number conservation in a MZ interferometer with phase gratings. 

A group in Colorado used Bragg diffraction at various orders to built MZ interferom¬ 
eter for a Ne* beam [99] (figure 12b). A Bragg scattering interferometer for Li atoms 
with high contrast and a count rate of 17 kc/s was used to measure the polarizability of 
Li atoms [101, 102, 103]. 


Talbot-Lau (near-field) interferometer. A high degree of spatial coherence is needed to 
create recurring self-images of a grating due to near-held diffraction (the Talbot effect). 
But completely incoherent light can still produce fringes downstream of a grating pair 
(the Lau effect). Two gratings with equal period d separated by a distance Li create a 
the Laue fringe with period d' at a distance L 2 beyond the second grating: 


(53) 


L2 


Li-L 


n 

2m 

n 


^ 2m 


d' = d 


L2 + Li 

mLi 


Here Lt = 2d?'/XdB is the Talbot length and the integers n and m refer to the revival 
of the order Fourier image. If a 3’’'^ grating with period d' is used as a mask to hlter 
these fringes, then a single large-area integrating detector can be used to monitor the 
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Fig. 13. - A sketch of the Talbot-Lau interferometer setup consisting of three gratings. The first 
grating is illuminated by an uncollimated molecular beam. Still, coherent interference occurs 
between all paths that originate from one point at the first grating and meet at the a point 
on the third grating. By varying the grating position 0:3, a periodic pattern in the molecular 
distribution can be detected. Figure from ref. [58]. 


fringes. In such a 3-grating Talbot-Lau Interferometer (TLI) the contrast is unaffected by 
the beam width and a large transverse momentum spread in the beam can be tolerated, 
hence much larger count rates can be obtained. The TLI does not separate the orders - 
components of the wave function are only displaced by one grating period at the Talbot 
length, it is still sensitive to inertial forces, decoherence, and field gradients. 

In a TLI the relationship Li = L 2 = means that the maximum grating period 

is d < 'jLiXdB ~ where M represents mass for a thermal beam. In comparison, 

for a MZI with resolved paths the requirement is d < XdBL/{Ax) ~ where Ax 

is the width of the beam and L is the spacing between gratings. Thus the TLI design 
is preferable for demonstration of interference for large mass particles (see lectures by 
M. Arndt). 

A Talbot-Lau interferometer was first built for atoms by John Clauser [104]. Using 
a slow beam of potassium atoms and gratings with a period of d =100 /rm, and a count 
rate of (!) = 4 x 10^ atoms/sec was achieved, even though the source brightness was 
2500 times weaker than in the 3 grating Mach Zehnder interferometer at MIT, but the 
signal was about 3000 times stronger. Because of its attractive transmission features, and 
the favorable scaling properties with XdB, the TLI has been used to observe interference 
fringes with complex molecules such as Ceo, C 70 , C 60 F 48 , and C 44 H 30 N 4 [60, 61]. 

2’2.2. Selected experiments with beam interferometers. 


Examples of phase shifts. We will now briefly discuss typical phase shifts that can be 
observed using an interferometer. In the JWKB approximation the phase shift induced 
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by an applied potential is given by: 

(54) 0(fco) = j {k{x) — fco(a:)) dx = j Ak{x) dx = ko j (n(x) — 1) dx 

Fci is the classical path and /cq and k are the unperturbed and perturbed fc-vectors, 
respectively and n is the refractive index as given in equation 31. If the potential V is 
much smaller than the energy of the atom E (as is the case for most of the work described 
here) the phase shift can be expanded to first order in V/E. If V is time-independent, 
one can furthermore transform the integral over the path into one over time by using 
t = xjv. 


(55) Aip(ko) = —— [ V(x)dx = —— [ V(t)dt 

kv ri 

where n is the refractive index as given in Eq. 31. 

A constant scalar potential V applied over a length L^t results in a phase shift 
(j>{ko) = ■ The power of atom interferometry is that we can measure these 

phase shifts very precisely. A simple calculation shows that 1000 m/s Na atoms acquire a 
phase shift of 1 rad for a potential of only V = 6.6x 10“^^ eV in a 10 cm interaction region. 
Such an applied potential corresponds to a refractive index of |1 —n| = 2.7 x 10“^^. Note 
that positive V corresponds to a repulsive interaction that reduces k in the interaction 
region, giving rise to an index of refraction less that unity and a negative phase shift. 

Equation (55) further more shows that the phase shift associated with a constant 
potential depends inversely on velocity and is therefore dispersive (it depends linearly 
on the de Broglie wavelength). If, on the other hand, the potential has a linear velocity 
dependence, as in [/ = /r • i (v x E) for a magnetic dipole /x in an electric field E, the phase 
shift becomes independent of the velocity v [105, 106, 107, 108]. Similarly, a potential 
applied to all particles for the same length of time, rather than over a specific distance, 
will produce a velocity independent phase shift (/ = — F J)*'’ V{t) dt. The latter is related 
to the scalar Aharonov Bohm effect [109, 110, 111]. 

Electric polarizability . By inserting a metal foil between the two separated arms, as 
shown in figure 14, an a uniform electric field E can be applied to to one of the separated 
atomic beams, shifting its energy by the Stark potential U = —aE'^j^. The static scalar 
ground-state polarizability can then be determined from the phase shift, 4 >, of the 
interference pattern by 

(56) [IF) (2S„), 

where V is the voltage applied to one electrode in the interaction region, D is the distance 
between the electrode and the septum, v is the mean velocity of the atomic beam, and 
Leff is the effective length of the interaction region defined as (^) Leff = / E^dz. 
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Fig. 14. - Measurement of atomic polarizability, (a) Schematic of the interaction region installed 
behind the second grating, (b) Measured phase shifts vs. applied voltage. The two different 
signs of the phase shift stem from the voltage being applied on either the left (open circles) or 
the right (filled circles) side of the interaction region (arm of the interferometer). The fit is to 
a quadratic and the residuals are shown on the lower graph. Figure from ref. [97]. 


For an accurate determination of electric polarizability, the three factors in equa¬ 
tion (56) must each be determined precisely. They are (1) the phase shift as a function 
of applied voltage, (2) the geometry and fringing fields of the interaction region, and 
(3) the velocity of the atoms. In ref. [112] the uncertainty in each term was less than 
0.2%. This allowed to extract the static ground-state atomic polarizability of sodium to 
Oipoi = 24.11 X 10“^"^ cm^, with a fractional uncertainty of 0.35% [112]. Similar precision 
has been demonstrated for aue by the Toennies group [113] and a^i with a precision of 
0.66% by the Vigue group [102, 103]. These experiments offer an excellent test of atomic 
theory. 


Refractive index. A physical membrane separating the two paths allows to insert a gas 
into one path of the interfering wave. Atoms propagating through the gas are phase 
shifted and attenuated by the index 


(57) 




The phase shift due to the gas, 


(58) 


(j){N,z) = {2TTNkz/kcm)Re[f{kcm)], 


is proportional to the real part of the forward scattering amplitude, while the attenuation 
is related to the imaginary part. Attenuation is described by the total scattering cross 
section, and this is related to Im[/] by the optical theorem 


(59) 


Att 

o-tot = T —Ini[/(fcc 

^cm 
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Fig. 15. - (Left) Phase shift 0 as a function of gas density N for different gas samples. (Right) 
Phase shift vs Fringe amplitude. The fringe amplitude is proportional to Figure 

from ref. [97]. 


Measurements of phase shift as a function of gas density are shown in figure 15. 

The ratio of the real and imaginary parts of the forward scattering amplitude is a 
natural quantity to measure and compare with theory. This ratio, 


(60) 


p{k) 


m) 

ln[A{N)/A{0)] 


Re[/(fc)] 
Im[/(fc)] ■ 


where A is the fringe amplitude, gives orthogonal information to the previously studied 
total scattering cross section. In addition it is independent of the absolute pressure in 
the scattering region and therefore much better to measure. 

The motivation for studying the phase shift in collisions is to add information to long¬ 
standing problems such as inversion of the scattering problem to find the interatomic 
potential V{r)^ interpretation of other data that are sensitive to long-range interatomic 
potentials, and description of collective effects in a weakly interacting gas [114, 115, 116, 
117, 118, 119, 120, 121, 122]. These measurements of p are sensitive to the shape of the 
potential near the minimum, where the transition from the repulsive core to the Van der 
Waals potential is poorly understood. The measurements of p{k) also give information 
about the rate of increase of the interatomic potential V(r) for large r independently 
of the strength of V{r). The real part of / was inaccessible to measurement before the 
advent of separated beam atom interferometers. 


Measurement of the Coherenee length. The coherence length of an atomic beam is 
related to its momentum distribution Up by eq.26 to ic = -^ = ^' 27 ^^ The longitudinal 
coherence length ticoh limits the size of optical path difference between two arms of an 
interferometer before the contrast is reduced. 
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Fig. 16. - Measurement of the coherence length of a Na atomic beam with de Broglie wavelength 
of AdB = 170 pm in an atom interferometer. The phase shift was realized using the electric 
polarizability of the atoms, by applying an electric field. 


The coherence length can be measured directly in an atom interferometer. Applying 
a classical potential in one path results in a phase shift </> and simultaneously in a spatial 
shift of the wave function relative to the other path by Sx = ^^dB- This allows a 
direct measurement of the first order coherence function + Sx)) = 

O’* {x)'ip{x + Sx) dx. An example of a measurement for a Na atomic beam is shown in 
figure 16. 

2’3. Einsteins recoiling slit: a single photon as a coherent beamsplitter. - . Up to now 
when using ligh fields to manipulate atomic motion we were using classical light which 
is described by a classical electro magnetic wave. We now discuss the other extreme: An 
experiment where a single emitted photon is used as a beam splitter 

In spontaneous emission an atom in an excited state undergoes a transition to the 
ground state and emits a single photon. Associated with the emission is a change of the 
atomic momentum due to photon recoil [123] in the direction opposite to the photon 
emission. The observation of the emitted photon direction implies the knowledge of the 
atomic momentum resulting from the photon-atom entanglement. 

If the spontaneous emission happens very close to a mirror the detection of the photon 
does not necessarily reveal if it has reached the observer directly or via the mirror. For 
the special case of spontaneous emission perpendicular to the mirror surface the two 
emission paths are in principle in-distinguishable for atom-mirror distances d c/T 
with c the speed of light and T the natural line-width. In this case the photon detection 
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Fig. 17. - Motional coherence generated by a single spontaneons emission event, (a) An atom in 
front of a mirror spontaneously emits a single photon. For emission perpendicular to the mirror 
surface an observer can in principle not distinguish if the photon has been reflected or not. 
Momentum conservation in the atom-photon system implies that the atom after the emission 
is in a coherent superposition of two different momentum states separated by twice the photon 
recoil, (b) With the spatial extension of the atom corresponding to the optical absorption cross 
section, indistinguishability can be estimated by the projected overlap of atom and its mirror- 
image. This overlap is represented color coded on a sphere for all emission directions (red: full 
coherence, blue: no coherence) and on the pattern at the atom detector, (c) For large distance 
between the atom and the mirror the coherent fraction drastically reduces, approaching the 
limit of vanishing coherence in free space. Figure from ref. [124]. 


projects the emission in a coherent superposition of two directions and the atom after this 
emission event is in a superposition of two motional states. Consequently the photon 
can be regarded as the ultimate lightweight beamsplitter for an atomic matter wave 
(figure 17a). Consequently spontaneous emission is an ideal model system to implement 
the original recoiling slit Gedanken experiment by Einstein [125]. 

This reasoning can easily be generalized to the case of tilted emission close to the 
mirror surface. One finds residual coherence for emission angles where the optical ab¬ 
sorption cross section of the atom and the mirror-atom observed by a fictitious observer 
in the emission direction still overlap (figure 17b). For larger distance to the mirror, the 
portion of coherent atomic momentum is strongly reduced (figure 17c). 

The coherence can be probed by superposing the two outgoing momentum states 
using Bragg scattering at a far detuned standing light wave on a second mirror [82, 
69]. One observes an interference pattern as function of a phase shift applied by 
translating the Bragg standing light wave by moving the retro-reflecting mirror. The 
two outermost momentum states, which represent maximum momentum transfere due 
to photon emission in the back to back directions orthogonal to the mirror surface are 
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Mirror position L [pm] 

2,58 2.75 2.92 3.09 3.26 




Mirror position L [pm] 



Fig. 18. - Experimental confirmation of coherence induced by spontaneous emission employed as 
the first beamsplitter of an atom interferometer. The recombination is accomplished by Bragg 
scattering from a standing light wave. The relative phase of the two paths can be changed by 
moving the ’’Bragg” mirror as indicated, (top) For a large mean distance between atoms and 
the mirror of 54 p.m no interference signal is observed confirming the free space limit. The inset 
depicts the position of the mirror relative to the atomic beam, (bottom ) For a mean distance 
of 2.8 pbm the two complementary outputs of the interferometer reveal an interference pattern 
with a maximal visibility of 5.9% ± 1.1%. Figure adapted from [124]. 


expected to show the highest coherence. 

The experiment of ref. [124] was performed with a well collimated and localized beam 
of ^°Ar atoms in the metastable IS 5 state (for the level scheme see figure ??) (figure 4). 
In order to ensure the emission of only a single photon we induce a transition IS 5 -T 2 p 4 
(A_e = 715 nm). From the excited state 2 p 4 the atom predominantly decays to the 
metastable IS 3 state via spontaneous emission of a single photon (Xse = 795 nm, branch¬ 
ing ratio of IS 5 /IS 3 = 1/30). The residual IS 5 are quenched to an undetectable ground 
state with an additional laser. Choosing the appropriate polarization of the excitation 
laser the atomic dipole moment is aligned within the mirror plane. The interferome¬ 
ter is realized with a far detuned standing light wave on a second mirror. Finally the 
momentum distribution is detected by a spatially resolved multi channel plate approx¬ 
imately Im behind the spontaneous emission enabling to distinguish between different 
momenta [124]. 

Figure 18 shows the read-out of the interference pattern for two distances between 
the atom and the mirror surface. The upper graph depicts the results obtained for 
large distances (> 54 ^m) i.e. an atom in free space. In this case no interference is 
observed, and thus spontaneous emission induces a fully incoherent modification of the 
atomic motion. For a mean distance of 2.8 ^m clear interference fringes are observed 













Interferometry with Atoms 


41 


demonstrating that a single spontaneous emission event close to a mirror leads to a 
coherent superposition of outgoing momentum states. 

It is interesting to relate this experiment to the work by Bertet et al. [126] where 
photons from transitions between internal states are emitted into a high fines cavity, 
there the transition happens from from indistinguishability when emission is into a large 
classical field to distinguishability and destruction of coherence between the internal 
atomic states when emission is into the vacuum state of the cavity. Using the same 
photon for both beamsplitters in an internal state interferometer sequence, coherence 
can be obtained even in the empty cavity limit. In this experiment the photon leaves 
the apparatus and one observes coherence only when the photon cannot carry away 
which-path information. This implies that the generated coherence in motional states 
is robust and lasts. In this sense it is an extension of Einstein’s famous recoiling slit 
Gedanken experiment [125]. In free space the momentum of the emitted photon allows 
to measure the path of the atom. This corresponds to a well defined motional state of 
the beamsplitter i.e. no coherence. Close to the mirror the reflection renders some paths 
indistinguishable realizing a coherent superposition of the beamsplitter. The large mass 
of the mirror ensures that even in principle the photon recoil cannot be seen. Thus the 
atom is in a coherent superposition of the two paths. 
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3. — Interferometry with Bose-Einstein condensates in double-well potentials 

It was recognized as early as 1986 by J. Javanainen [127], almost 10 years before 
the experimental production of the first Bose-Einstein condensates (BECs) in dilute 
gases [128, 129], that BECs trapped in double-well potentials shared common features 
with solid states Josephson junctions [130]. In the former case, the cooper pairs are 
replaced by neutral atoms, and the thin insulating layer by a tunnel potential barrier 
that can be realized either optically [131, 132, 133, 134], magnetically [135], or with 
hybrid traps such as radio-frequency dressed magnetic potentials [136, 137, 138, 139, 140]. 
Nevertheless, the contact interaction between atoms, and the absence of leads connecting 
to an external circuit modifies the physics compared to standard Josephson junctions. 

In section 3'1, we start by neglecting interactions and introduce the basic concepts 
associated to this system: its reduction to a two-level problem (3" 1.2), its dynamics 
(3’1.3), and the ways to measure it and use it for interferometry (sections 3'1.5, 3'1.6 
and 3'1.7). In section 3'2, we analyze the effects arising from interactions: the emergence 
of the nonlinear Josephson dynamics (3'2.2), the possibility to control the quantum 
fluctuations (3’2.3) by splitting a condensate (3'2.4) and connections to interferometry. 

3’1. A Bose-Einstein condensate in a double-well potential: a simple model. - Let 
us first assume that interactions between particles are negligible, and that a condensate 
containing N atoms is trapped in a one-dimensional (ID) tunable potential, which can be 
turned into a double-well potential [141]. We assume that the temperature is negligible, 
i.e. that the system is initially in the ground state. In practice, experiments generally in¬ 
volve interacting atoms at Hnite temperature, which are trapped in real three-dimensional 
potentials, but the above simplifications will help us tackle the most important features 
of the true system. The effect of interactions will be the subject of section 3’2. The 
geometry is depicted in figure 19. 



Fig. 19. - A simple model describing a BEC in a double well potential. The horizontal lines 
correspond to the single-particle eigenenergies. The two lower lying orbitals are also depicted, a, 
In the single-well limit, the states are approximately equally spaced in energy, b, In the presence 
of two minima, the lower-lying states gather by pairs, which become completely degenerate in 
the strongly split limit, c, A slight asymmetry of the trap is enough to lifts the degeneracies, 
and the two lower-lying states are localized in the wells. 


































Interferometry with Atoms 


43 


3’1.1. Single-particle approach. Upon neglecting interactions and temperature effects, 
the condensate wave function (the ground state of the A^-particle system) can be written 
as a product of single particle ground states: 


N 

(61) ^'BEc(a;i,a^2, • ■ ■,xn) = 

i=i 


Here, 'tpg is the single particle ground state, i.e. the solution of h-tp = Etp having the 
smallest possible energy E = Eg. The single particle Hamiltonian h is given by 


(62) 


h 


2 m dx^ 


V{x). 


In a word, because the system is non-interacting and prepared at zero temperature, it is 
formally equivalent to N independent single particles prepared in the same initial ground 
state of the double well potential. We can thus understand everything by studying the 
case of a single particle trapped by the potential of figure 19. 

The spectrum of h is displayed in figure 19. When the two wells are not separated 
and the barrier height small, the spectrum resembles that of a harmonic oscillator, that 
is, the levels are roughly equidistant with a splitting E ~ hui, where w is typically the 
trap frequency. If one is able to prepare the system at a temperature T <C E/Rb, the 
system is almost purely in the ground state. When the spacing between the traps (or 
the barrier height) is increased, the eigenstates group by pairs, until the pairs become 
completely degenerate which corresponds to a probability to tunnel from one side to the 
other that has become negligible (cf. Hgure 19). 

3T.2. Two-mode approximation. Guided by the fact that the system can be prepared 
in the ground state, we here simplify the full problem (Schrodinger equation with the 
Hamiltonian (62)) to an effective one, which will describe the physics at low energy. 
The most simple description taking into account deviations from the ground state is 
obtained by keeping only the two lower lying states ipg and i/'e- In the next sections, we 
will distinguish different cases depending on the precise shape of the trapping potential 
V{x). We will first consider a strictly symmetric trap: V{—x) = V{x), and then discuss 
the differences arising from an asymmetry. We note also that one can go beyond the 
two-mode approximation [142], but this goes beyond the scope of these lectures. 


Symmetric traps. The symmetry with respect to the barrier ensure that each eigenstate 
of the system is itself symmetric or antisymmetric. Eigenstates thus correspond to an 
equal probability of Hnding the particle in one well or the other (the particles are not 
localized around one potential minimum as one would expect in classical physics). It 
seems more natural to describe the system in terms of modes that are localized on each 
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side of the barrier, consequently we introduce the left/right basis(^^): 

(63) ^ ^ ’ 

(64) il,^= ^ - ij ,^). 

This corresponds to a 45° rotation in the 2-dimensional Hilbert space spanned by the 
initial basis {■0g,^e} describing the problem. 

From the Hamiltonian (62), we can derive a simple effective Hamiltonian describing 
the system in the left/right basis introduced above. We approximate the wave function 
1 /( 0 ;, t) by the Anstaz 

(65) -ipix, t) = + Crit)'lprix) 

which corresponds to a two-mode approximation. Following the approach of refs. [127, 
143, 144, 145] we obtain the two-mode Hamiltonian 

( 66 ) "/) 

for the evolution of the state vector ^ = (ci, Cp). The matrix elements of this Hamiltonian 
are 

(67) E = (f/ll^l^l) = {'lpr\h\'tpr) = ]^{Eg + Ee), 

(68) J = -{HMA) = -{Ami) = liEe- Ag), 

Eg^e being the eigenenergies of the ground and first excited eigenstates respectively. 

Sensitivity to trap asymmetry. A problem arising when splitting a non-interacting gas is 
the sensitivity to the asymmetry of the trap. If we assume that the right well is higher 
by an energy i5 > 0 compared to the left well, the Hamiltonian (66) becomes(^^) 

(09) '-=(4'/ f/)^ 


(^^) We have chosen the wave functions V’g)®) S'lid ipeix) real, and with the proper sign, such 
that i/>i(a;) and '4>-r{x) as defined by equations (63) and (64) are indeed localized around the left 
and right wells respectively. 

(^^) Compared to equation (69), we have removed the constant energy shift E which does not 
play any role. 
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In the absence of asymmetry (<1 = 0), the eigenstates are the symmetric and antisym¬ 
metric superposition of left and right modes as explained above. Nevertheless, a small 
asymmetry is responsible for the collapse of the eigenstates in the two wells. This situ¬ 
ation was studied in a box potential in reference [146], but the discussion below will be 
valid for any potential. The eigenvalues of (69) are 

(70) = 


where the “-I-” corresponds to the excited state, and the ” to the ground state. Here 
we have introduce an auxiliary angle 9 defined by tan0 = —2JI6. The eigenstates are 
given by (see ref. [147], Volume 1, Complement Biy) 


(71) 

(72) 


V’g = cos -ipi + sin 

f/'e = - sm -Ipx + cos -ip^. 


We see that if there is no coupling, J = 0, and in the absence of a bias, 5 7 ^ 0, then 9 = 0 
and the eigenstates ipg and i/ig “collapse” to the left and right wells respectively. From the 
expression of the eigenvectors, one could also compute the population of the two modes 
as a function of the ratio J/S. This sensitivity of the double well to an asymmetry is 
illustrated in figure 19c, where a small gradient has been added to a symmetric double¬ 
well potential, leading the two lower lying states to collapse into the left and right wells. 

3’1.3. Rabi dynamics. We now focus on the dynamics that can occur when the system 
is not prepared in an eigenstate, or when the potential is changed such that the eigenstates 
are modified, which brings the system out of equilibrium and triggers dynamics. In the 
left/right basis, we recognise in ( 66 ) a simple two-level Hamiltonian with off-diagonal 
coupling, which will exhibit oscillations of the populations |ci|^ at the (angular) Rabi 
frequency 


(73) 


Qr = ijjh. 


If the system is initially prepared in a superposition of the two eigenstates, the par¬ 
ticles will tunnel back and forth between the two wells at this frequency. The corre¬ 
sponding time tr = h/2J is the typical time scale associated to tunnelling of atoms. 
This oscillation of the populations is accompanied by an oscillation of the relative phase 
(j) = arg(cj'cr) at the same frequency. While the oscillation of the populations is strictly 
sinusoidal: |ci(t)|^ = [l-|- 74 cos(ORt-|-(p)]/ 2 , with —1 < H < 1 , the oscillation of the phase 
becomes anharmonic as soon as the amplitude \A\ is not much less than 1. Note that 
when the initial state is localized in one of the wells, the amplitude is maximal [Hj = 1 , 
which means that all the atoms tunnel back and forth between the two wells. 


3’1.4. Splitting. 



46 


J.-F. SCHAFF, T. Langen and J. Schmiedmayer 


Principle of coherent splitting. From the results of section 3'1, we see that a conden¬ 
sates can be coherently split by adiabatically changing the potential from a single well 
to a double well. This was for example studied numerically in ref. [148] and realized 
experimentally [149, 136, 139, 150, 140]. Adiabaticity ensures that the system, initially 
prepared in the ground state, will stay in the ground state during the transformation, 
which, analysed in the right left basis, is a coherent superposition state of the type 

(74) -tf;= + . 

In order for the final state to be a coherent superposition after splitting, the additional 
constrain that the ground state must be delocalized in the two wells must be fulfilled. 
As explained in section 3T.2, this condition is quite constraining, as even a very slight 
asymmetry would lead the ground state to collapse in the well having the lowest energy. 
This is modified in the case of an interacting Bose-Einstein condensate, in which a new 
energy scale related to interactions, the chemical potential, will play a significant role. 


Condition for adiabatic splitting. We model the splitting process as a slow and monotonic 
change of the parameters J{t) and 6 (t) which are now time-dependent. The condition 
for adiabatic evolution is given by [151] 


(75) 


hillJeldtllpg) 

Eg — Eg 


Here the two eigenstates of (69) have the energies given by equation (70). When evalu¬ 
ated, condition (75) yields 


(76) 


jj + S/2\ 
4J2 -t,52 


< 1 . 


In the symmetric case discussed above {S = 0) we see that, during splitting, because 
J —> 0, the splitting must become slower and slower in order to stay adiabatic. It will 
not be strictly speaking possible to stay adiabatic all the way to J = 0. 


Non-adiabatic splitting. This opens the possibility of splitting non-adiabatically: if the 
ground state wave function is still delocalized when adiabaticity starts to break down, 
then the wave function will stay essentially the same as the splitting process goes on 
despite the fact that the eigenstates may change very abruptly. In this case, a coherent 
superposition of left and right modes is also created. 

3T.5. Time-of-flight recombination. A coherent superposition state can be measured 
by recombining the two “arms” of the interferometer (in our case the two BECs trapped 
in the two wells) and measuring the interference. Time of flight recombination is similar 
to a double split experiment. The matter-waves in each potential minima are released and 
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Fig. 20. - Principle of time-of-flight recombination. After a long time of flight, the two wave 
packets overlap and interfere. Figure adapted from Tim Langen’s PhD thesis [152]. 


freely expand until they overlap, forming a typical double-slit interference pattern (cf. 
examples of such patterns on figure 21). Nevertheless, as surprising as it may seem, the 
observation of interference fringes between two overlapping Bose-Einstein condensates 
is not a proof of coherence. Indeed, it was shown that two condensates having each 
a fixed, well defined, atom number (Fock states), and therefore no fixed relative phase 
(vanishing off-diagonal terms in the one-body density matrix) would display interference 
when overlapped and measured [153]. The fringes that appear can be seen as the result of 
the measurement-induced collapse of the system’s state to coherent state (cf. section 3'2) 
having a random relative phase. Therefore, first order coherence is proven when the phase 
can be measured repeatedly and always yields the same value, that is when the fringes 
observed are identical from shot to shot. 

In the current and next sections (sections 3T.5 and 3T.6) we assume that a condensate 
has been split symmetrically, and that a relative phase (j) has been accumulated. The 
initial wave function can therefore be written 

(77) t = 0) = + e*‘^V’r(a;)) 

with 'ipT-{—x) = '0i(ai) = f{x,t = 0)(^^). The generalization to an imbalanced splitting 
yielding different populations in the left and right modes is straightforward. 

When the two trap minima are well separated, one can generally approximate them by 
harmonic potentials, and the left and right modes are therefore well modeled by minimal 
uncertainly Gaussian wave packets. Let us assume a symmetric trap with trapping 


(^^) Note that in quantum mechanics, the global phase of the state (77) can be chosen arbitrarily. 
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frequencies w. When the trapping potential is abruptly turned off, the dynamics can be 
decomposed into 

1 . the center of mass movement, which will follow a classical trajectory: Zcom = 
-9tV2, 

2 . the relative movement in the free-falling frame, witch is identical to that in the 
absence of any gravity. 

We can thus neglect gravity. For simplicity, we assume that the 3D potential is separable, 
and just solve the problem in the splitting direction x. 


Expansion of a single wave packet. The function / is initially 

with a;i = d being the distance between the trap minima. Equivalently, in mo¬ 

mentum space representation(^‘^) 

(79) f(p.t = 0) = exp exp {-if) . 

Here (Jq = \/h/mui is the spatial width of the ground state in each well. The width in 
momentum space is h/a^ = \/mhuj. The free expansion after trap switch off is described 
by the evolution operator tj{t) = exp{—iHt/h) = exp{—ip^t/{2hm)), which is diagonal in 
momentum representation thanks to the absence of a potential in H. The wave function 
is therefore given at time t by multiplying the initial one (79) by U and applying an 
inverse Fourier transform to recover the spatial wave function. We recover the well 
known results that a Gaussian minimal uncertainty state stays Gaussian with a width 
increasing as 

(80) a{f) = (Jo (l -I- . 


The time-dependent wave function is given by 


(81) 


f{x,t) = 


1 


1/4 


^ \ 2a{tY 

with the spatially-dependent phase accounting for the velocity distribution of the atoms 


{x - Xif 


exp{iip{xf)), 


(82) 


(p{x,t) = 


ht 


2to(Tq 


{x - Xlf 

aitr 


- 1 


(!'*) We have chosen the definition of the Fourier transform and its inverse as g{p,t) = 
f dx, and g{x,t) = (27rh)~^''^ f t) dp. 
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Fig. 21. - Collection of Bose-Einstein condensates interference patterns, a, First BEC interfer¬ 
ence pattern observed by Andrews et al. [149]: the BEC was trapped magnetically and split 
with a repulsive optical barrier, b, c, d, e interference pattern of Schumm et al. [136], Jo et 
al. [139], Hofferberth et al. [150] and Berrada et al. [140]: the BECs were trapped and split in 
radio-frequency dressed magnetic traps. The patterns of c and d display a relative phase which 
varies along the long axis of the BEC (see section ?? for a discussion of such ID effects). 


Interference of two expanding wave packets. In the case of a double-well potential, the 
two wave packets follow the evolution described above, and, using (77) as the initial 
state, we obtain an interference pattern when they start to overlap. The interference 
term is 


(83) 


t)\\f{-x, f )I cos {2k{f)x -h (/>), 


where the time-dependent wave vector k{t) is 


(84) 


htd 

2mcrg(T(t)2 


At long times {cat ^ y/(4 (Tq — 1), the width of the wave functions is much larger than 
their initial distance d, and the two envelopes coincide f{—x,t) ~ f{x,t). The density is 
approximately 


(85) 


N\tp{x, t)p ~ iV|/(a;, f)p (1 -|- cos {2k{t)x + (f)). 
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Position 


Fig. 22. - Effect of mean-field repulsive interaction on the fringe spacing: The atom wave 
moving expanding from one well moves throngh the expanding atom wave of the other well and 
accumulates a phase shift. Consequently the fringe spacing becomes position dependent, leading 
to bent fringes. 


We see that the fringe spacing Ap = ■n/k{t) coincides with Ap = AdB/2, half the de 
Broglie wavelength of an atom initially position at a trap minimnm, and having a velocity 
V = ±d/2t such that it could reach the central point a: = 0 after a time t. 

3’1.6. In-trap recombination. Reading out the relative phase of two waves by overlap¬ 
ping them and observing the interference pattern, as described above, is not the standard 
way used in optics. The beams are rather recombined on a beam splitter. The two input 
modes are “mixed” and the resulting two output modes have different intensities, depend¬ 
ing on the relative phase of the input modes. A beam splitter is therefore an object which 
converts phase differences into intensity differences and vice versa. This strategy can by 
extended to Bose-Einstein condensates trapped in double-well potentials [148, 140]. 

Recombining with Rabi oscillations. A first possibility is to adiabatically reconple the 
two traps by reducing the barrier height or well spacing [148]. In this case the cou¬ 
pling J, which is initially negligible, increases, and Rabi oscillations are triggered (cf. 
section 3T.3). A quarter of a Rabi oscillation will convert the phase into a population 
imbalance between the two wells, which can then be measured. This requires to couple 
the wells for a time r ~ h/8J (assuming that Rabi oscillations have not started during 
the recombination stage). 

Nevertheless, we will see that interactions will reduce the efhciency of this approach, 
as they tend to reduce the amplitude of the oscillations of the populations (see section 
3’2.2, equation (108)). 

Non-adiabatic recombination. An alternative approach is to modify the trap in a non- 
adiabatic way. Here we rely on the symmetry of the trapping potential V{—x,t) = 
V{x,t). Before recombination, the wave function describing the two condensates has the 
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Fig. 23. - Accumulation of a relative phase due to an energy difference 5 yf 0 (cf. equation (89)). 
The phase is measured by time of flight recombination. Inset: evolution of the interference 
pattern (see figure 21e) from which the phase is measured. Figure adapted from ref. [140]. 


same symmetry, except for a relative phase (cf. equation (77)). From the symmetry of 
the potential and propagating with the Schrodinger equation, one finds that the BEC 
wave function can be written at all time as 




( 86 ) 


The imbalance 2 = |'0P dx — dx can readily be calculated from Eq. (86), 

which gives 


(87) 


z{t) = C{t) sin^. 


Note here that the argument is very general, and has not assumed a particular time- 
dependence or shape of the potential V{x,t). It only relies on the symmetry of the 
potential and wave function, and on the superposition principle. The contrast 



( 88 ) 


depends on the details and can be maximized by optimizing both the shape and the 
way the potential is dynamically varied. As a whole, this sequence turns an initial state 
which has equal population on both sides of the barrier and a relative phase 4>, into a 
state whose populations depend on the sine of the phase 4>. This is very similar to an 
optical beam splitter. 

Note that other strategies have been attempted, such as measuring the phase-dependent 
heating produced when two condensates are merged [154]. Nevertheless, this technique 
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Fig. 24. - Mach-Zehnder interferometer with a Bose-Einstein condensate in a double-well poten¬ 
tial. a, Nearly adiabatic splitting creates a coherent superposition of left and right modes. The 
recombination is performed non-adiabatically. b, At the end of the sequence, the atom number 
difference between the two traps show interference fringes. Grey dots: single realizations, black 
dots: ensemble average. Figure adapted from ref. [140]. 


strongly relies on interactions, as the heating was interpreted in terms of soliton creation, 
which could not occur in the absence of interactions. 

3’1.7. Phase shifts. Here we assume that the two wells are not coupled (J = 0), 
and that the condensate has been split symmetrically. If a force F is applied along 
the splitting direction during a time t (or equivalently if the system is subject to an 
acceleration a = —F/m), the energy shift <5 = Fd between the two wells gives rise to a 
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phase shift, given at first order by 


(89) 


Fdt 

^force — Z 5 

h 


where d is the distance between the two trap minima. This illustrates how such an atom 
interferometer could in principle be used to measure inertial effects. An example of such 
phase accumulation is displayed in figure 23 for different values of 5. This assumes that 
both traps stay harmonic with the same frequencies. If they have a different trapping 
frequency wi = Wr + Aw, then there is an additional phase shift 


(90) 


^asym. — 2^^^' 


This illustrates the difficulty of performing precision measurements with such interfer¬ 
ometers, because this would require a very good knowledge of the trap shape to know its 
response to an external force. 

3’2. Effect of interactions: Josephson dynamics, squeezing, and dephasing. - In this 
section, we extend the results of section 3T to the case of interacting particles. A simple 
way to take into account interactions is to still assume that the particles are distributed 
among two modes, say the left and right modes introduced in section 3T.2, and to work 
with the second quantization formalism. To each mode is associated a pair of creation 
and annihilation operators. It is convenient to work in the basis of Fock states, in which 
the number of atoms in each mode is well defined (i.e. has no quantum fluctuations). Let 
us write |ni,ni.) the Fock state having ni atoms in the left-localised mode and rir atoms 
in the right-localized mode. For example, the vacuum state (no atoms) is written |0,0). 
If we assume that the total atom number N is fixed, then the Hilbert space has -|- 1 
dimensions. This can be seen because the set of Fock states {\N — k,k),k € {0, A}} 
contains A -|- 1 vectors, and is an obvious basis. If we write I and P the annihilation and 
creation operators associated to the left mode, and f, those associated to the right 
mode, we have the usual relations 


(91) 

(92) 

(93) 

(94) 

(95) 


[i,P] = 1, [r,fp = 1, =0 

l\ni,n„) = ffrTilni - l,nr) 

P\ni , nr) = Vni + l|ni -I- 1, n^) 
r\ni,n,) = - 1) 

rpni,nr) = Vn7+T|ni,nr -|- 1) 


and i|0,nr) = 0, f|ni,0) = 0. The operators Pi = hi and fff = ff give the number of 
atoms in the left and right modes respectively. 
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3’2.1. Bose-Hubbard model with two sites. Once again we focus on the most simple 
model which will display the basic features arising from interactions: the Bose-Hubbard 
model. 

It has become increasingly popular in the past decade, thanks to the development of 
optical lattices in cold atom experiments. A double well potential can be seen as the most 
simple lattice: one containing only two sites. The usual lattice approach will therefore 
breakdown, due to the very small size of the system, but we will see that the basic 
features known in the case of large lattices, such as the existence of a superfluid phase 
and a Mott insulator phase, have an analogy in the case of only two sites. To properly 
derive the Bose-Hubbard model as an effective theory describing the real system (3D 
trapping potential), one can use a two-mode approach as in section 3'1, but this time on 
the field operator describing the system of interacting bosons. The field operator is thus 
approximated by 

(96) i’ix) = i'ipiix) + riprix)- 

We refer the interested reader to refs. [143, 144, 145, 155] for an example of such a 
derivation. 


Bose-Hubbard Hamiltonian. The Bose-Hubbard Hamiltonian reads 


(97) 


Hbh = — J 






hi(hi - 1)-I-hr(hr - 1) +^inr-ni) 


where the new energy scale U, which did not appear in section 3T, is related to inter¬ 
actions. The first term is a kinetic energy, and allows tunneling between the two wells. 
The third term corresponds to an energy difference S between the two modes (in the 
absence of tunneling), and the second term corresponds to the interaction energy associ¬ 
ated with contact interaction between the atoms (s-wave scattering). We note that this 
interaction is local, i.e. the atoms in the left-localized well do not interact with those in 
the right-localized well. In the following, we will always assume repulsive interactions 
[/ > 0 . 

In order to relate U to a physical value that can be measured in the experiment, we 
can analyze the Hamiltonian (97) in the absence of tunneling and bias (J = S = 0). 
The energy of the condensate in the left well is then given by Ei = Un\{n\ — l)/2. The 
chemical potential of this BEC is thus 

(98) fii = ^ ~ Um, 


therefore we see that the Bose-Hubbard model assumes that the chemical potential scales 
linearly with the atom number, and we can also make the identihcation 


(99) 


dN 


N—ni 
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being the chemical potential of the left-localized condensate, seen as a function 
of its atom number(^®). Therefore an underlying assumption in equation (97) is that 
the trap is sufficiently symmetric and the populations not too different, such that the 
interaction energy U is the same in both wells. 

3’2.2. Mean-field treatment. To get a first intuition on this Hamiltonian, we can 
apply a mean field approach to obtain the approximate dynamics of mean values of 
the operators [143, 144, 145, 155]. For that, one replaces the annihilation operators by 
complex numbers: 


(100) ( —)■ Cl € C, r —>■ Cr G C 

and the creation operators by their conjugate: cf, —)■ c*. The two c-numbers ci 

and Cr play the role of a (single particle) wave function, which is just defined on the two 
sites. This wave function is normalized to N: \c\\^ + jcrp = N, such that jcip ~ (hi) = ni 
and jcrp ~ (hr). We can write these numbers in polar representation: 


( 101 ) 

( 102 ) 


Cl = 7hle*'^‘, 

Cj. = 


It is then convenient to define the population imbalance z and relative phase (j) of the 
Josephson junction as 

(103) z = ^{ni-n,) = c*Cr), 

(104) (j) = (j)^ - (1)1 = arg(c(Cr). 

Substituting equations (100), (101) and (102) in the Hamiltonian (97), and rewriting it 
in terms of the variables introduced above, yields the classical Hamiltonian)^®) 

(105) Hmf = “ IIrVi - cos(/). 

n 2h 

In the absence of bias 5 = 0, it resembles the Hamiltonian of a pendulum with variable 
length in which z plays the role of momentum, and (j) is the angle with respect to the 
acceleration of gravity g: H = /{2m) — mg£ cos (j), i being the length of the pendulum, 
and m its mass. A difference is that the length £ gets shorter with higher momenta: Hmf 
is the so-called “momentum-shortened” pendulum. 


(^®) Note that the expression (98) cannot exclude the expression U = pi/ni compared to equa¬ 
tion (99). An answer to this ambiguity is obtained by using the full many-body Hamiltonian of 
a EEC in a 3D potential. In this case, the chemical potential does not scale linearly with m, 
which allows to check that equation (99) is the right form. 

(^®) We have use the relations: (hi) = (1 -|- z)N/2, (hr) = (1 — z)N/2, (J (hi)(hr) = Vl — z^N/2, 
and A = hi -I- hr. 
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Fig. 25. - Phase portraits of the bosonic Josephson junction illustrating Rabi oscillations (7 = 0), 
Josephson oscillations (7 > 1/N, blue lines) and macroscopic self-trapping (7 > 2/N, red lines). 
The red dots indicate the fixed points (see text). Figure adapted from T. Berrada’s PhD 
thesis [156]. 


The evolution of z and </) is obtained through Hamilton’s equations (/> = dHyiF/dz, 
and i = —dHuF/dcj), that is 


(106) 


z = —Hra/I ~ 


■ Hr 


Vi—2 


: cos (j). 


The dynamics described by these coupled nonlinear equations has been extensively stud¬ 
ied theoretically [143, 144, 145, 157], observed in experiments with BECs trapped in 
double-well potentials [132, 133, 140], and beautifully verified with BECs in two different 
hyperfine states [158] (^^). 

We shall just underline the existence of three distinct regimes, depending on the 
value of 7 = U/2J. These regimes are apparent of figure 25, which shows the phase 
space trajectories described by the system of equations (106)(^®): 

1. for 7 <C 1/N, interactions are essentially negligible, and we recover the Rabi dy¬ 
namics described in section 3T.3. In particular the angular frequency of these 
oscillations is the Rabi frequency Hr = 2J/ti. Close to 7 = 1/N, the oscilla¬ 
tion frequency around the fixed point {(j) = tt, z = 0) is significantly modified: 
w~ = HrVI - 'jN. 


(^^) In this experiment, two BECs share a common spatial wave function, but the two modes 
correspond to two different hyperfine states that are “tunnel coupled” by a two-photon transi¬ 
tion. If the left and right modes are replaced by the two hyperfine states, the system is formally 
equivalent, as it is described by the same Hamiltonian. The tunnel coupling energy J corre¬ 
sponds to the Rabi frequency of the two-photon transition, and interactions can be tuned thanks 
to a Feshbach resonance. 

(^®) For simplicity here we set the bias to <5 = 0. 
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2. for 1/N < 7 < 2/N, the fixed point (^ = tt, z = 0) becomes unstable, and two 

new fixed points appear at {(j> = Tr,z = ±a/ 1 — Depending on the initial 

conditions (^( 0 ),z( 0 )), the population imbalance and phase will either oscillate 
around 0 (so-called “Josephson oscillations”, blue trajectories in figure 25) at the 
“Josephson” or “plasma” frequency 

(107) ccp = DrV1 + 71V, 

which is higher than the Rabi frequency, or around the two other fixed points. In 
this later case the oscillation frequency is — (yfV)”^. 

3. for 7 > 2/N the modes discuss above still exist, but a new behavior appears, in 
which the phase is not a periodic function of time (it increases almost linearly), and 
the population imbalance oscillates around an non-zero value (red lines in figure 25). 
This is the so-called “macroscopic quantum self-trapping” regime. Interactions 
inhibit full population transfer. There is a smooth transition to the trivial case of 
two decoupled coherent condensates (7 —>■ 00 ), for which the phase accumulation 
is strictly linear, and the population imbalance constant, as atoms cannot tunnel 
anymore. 

Maximal amplitude of the Josephson oscillations. The fact that a bifurcation separates 
Josephson oscillations from macroscopic quantum self-trapping also imposes a maximal 
amplitude to the Josephson oscillations, given by 

(108) Zmax = max \z{t)\ = ^ . 

t jjy 

We see that if Rabi/Josephson oscillations are used to perform a beam-splitter operation, 
as discussed in section 3T.6, then interactions will limit the contrast of the beam splitter 

to Zniax- 

3’2.3. Fluctuations and interferometry. 

Shot noise limit. Fluctuations play an important role in interferometry as they are re¬ 
lated to the fundamental sensitivity attainable. To illustrate that, let us assume that N 
photons are sent into an interferometer one after another. Let us assume that the mean 
number of photons that escape through one output, say the “left output port” L, is given 
by (^ 1 ) = (1^/2)(1 + sin^), (j) being an adjustable phase shift. If we set </> to 0, then on 
average, 50% of the photons escape in each output port. This means that the photons 
are each in an equal superposition of being in both outputs of the interferometer. If 
one tries to detect them, they will be found with a probability p = 1/2 in each output. 
For one photon, the variance of the number of photons which has escaped through L is 
Var(hp^) = p X 0^ J- p X 1^ — p^ = 1/4, the last term being the mean photon number 
squared. For N photons, the experiment is repeated independently N times, therefore 
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the variances just add up, such that Var(hi) = N x Var(h|^^) = iV/4 = Anf. The sen¬ 
sitivity to a phase shift, giving the ability to distinguish between nearby values of the 
phase, is 


(109) 


Ahi 1 

d{n\)/d(j) y/N' 


This is the so-called shot noise limit of interferometry, also called “standard quantum 
limit”, or just “standard limit”. Here we have derived it assuming a perfect interferometer 
(of maximal contrast), and independent particles(^^). We see that if one is now able to 
introduce correlations between the particles, one may be able to improve the sensitivity. 
For instance, if the probability of hnding a particle in the “right” output gets larger when 
a particle has escaped through the “left” one, the sensitivity may be brought below the 
shot noise limit: 5(j) < l/v)/V. 


Mapping to a spin N/2 system. In order to analyze the fluctuations properties of the 
Bose-Hubbard model, we first introduce the mapping of the system to an ensemble of 
N spin 1/2. For that we define new operators form the creation/annihilation operators 
introduced above: 

( 110 ) = + 

( 111 ) ^y = Yi ~ 

( 112 ) = = ^{hi-h,). 

One can readily check that they satisfy spin commutation relations: [ISPS'™] = zSp, 
where {l,m,p) is any circular permutation of {x,y,z). 


“Coherent” atomie states. Among the many quantum states that such a complex system 
can have, a class of states is particularly interesting: the coherent states. They can be 
parametrized by two angles, setting the mean direction in which the spin is pointing. 
They are defined as the eigenstates of a spin component in the (d, </>) direction, i.e. the 
eigenstates of = Sx sin 9 cos (j) + Sy sin d sin </> -|- Sz cos 9, with eigenvalue N/2, where 
9 and (j) denote the polar and azimuth angles [159]. They can be written 

(113) |6»,(/) =-^=^= ^cos^/^-be*'^sin^f1'^ |0,0), 


(^®) Note that other values of the phase would give worse sensitivities, as the fluctuations would 
be reduced, but the slope of (hi((^)) as well. 
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z 



AJ 


Fig. 26. - Quasiprobability of a number-squeezed state with squeezing factor ^jv ~ 0.2 rep¬ 
resented on the Bloch sphere. The z-axis corresponds to the number difference. Figure from 


ref. [160]. 


form which underlines the fact that these states are generated by “stacking” particles 
one after another into the same (single particle) coherent superposition. 

To better understand these states, we can calculate the number distribution corre¬ 
sponding to a coherent state pointing on the equator of the Bloch sphere. By developing 
the expression (113) we obtain the probability of finding k particles in one wells and 
N — k in the other 



(114) 


We recognize a binomial distribution of atoms in the two wells. 


Bloch sphere representation. The possible states of a single spin 1/2 are parametrized by 
the polar and azimuthal angles and can therefore be represented as points on the surface 
of the Bloch sphere. This representation can be extended to a spin N/2. For this, one 
can represent the many-body state I'F) using the quasiprobability distribution [159] 


(115) 


Q(0,<^) = |(0,<(.|vE-)P, 


which can be plotted on the surface of a “generalized” Bloch sphere. 
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Bose-Hubbard Hamiltonian. Rewriting the Bose-Hubbard Hamiltonian (97) with the op¬ 
erators (110)-(112) we obtain(^°) 

(116) /7bh = -2JS, + USl - 6S,. 

This from is particularly appealing, because it has a nice interpretation in terms of 
trajectories on the Bloch sphere. Indeed, if one considers the first term only (setting U = 
i5 = 0), it generates an evolution operator U(t) = exp{2iJSxt/h), that is a rotation around 
the X axis of the Bloch sphere at the angular frequency given by the Rabi frequency 
Hr = 2J/h. Similarly, the third terms drives rotations around the 2 axis at an angular 
frequency 5/h (phase accumulation). The interaction term (second term) is different, as 
it is cx S'^. Classically, it would be a rotation around the z axis, whose frequency depends 
on the position on the z axis (making the approximation S'f ~ (S'z)*S'z). Contrary to the 
other, this non-linear term, which stems from interactions, is responsible for non-trivial 
modification of the state and leads to dephasing, squeezing, and the generation of strongly 
entangled states [159, 161]. 

Fluctuations. For simplicity, let us consider the coherent state |7r/2, 0), i.e. the eigenstate 
of Sx. All the other coherent states can be transformed into this one by applying an 
appropriate rotation. By construction this state has a length N/2 along the x-axis and 
no fluctuations: /S.S1 = 0. Nevertheless, it has fluctuations along the other axes y, and 
z, that can be seen as a consequence of adding N 1/2 spins having the same direction 
(and strong fluctuations). The uncertainty relations associated with the commutation 
relations satisfied by the components of the spin are ASk^Sm > l('§n)|/2. Applied to 
the coherent state above with {k, m, n) = {y, z, x) gives ASy ASz > N/A. Coherent states 
are characterized by a saturation of the uncertainty relation above and equal fluctuations 
along these two orthogonal axis: ASy = ^/N /2 and ASz = /2. Since Sz corresponds 

to half the population difference (cf. equation (112)), we see that coherent states have 
binomial number fluctuations Ah = If we express the fluctuations along y as 

fluctuations of the angle (j>, we obtain A(f> ~ {VN/2)/{N/2) = l/'/N. This is assuming 
Acj) <C TT, i.e. that the atom number be not too small. 


Squeezing as a resouree for sub-shot-noise interferometry. Squeezing is defined as a redis¬ 
tribution of the fluctuations along two axes orthogonal to the mean spin direction [159]. 
One can thus talk about number squeezing when the number fluctuations are reduced 
compared to the binomial case. For a state close to the coherent state l7r/2, 0)(^^) this 


(^®) We have assumed a closed system with constant atom number N = h\ -\- hr, and removed 
the constant terms from the Hamiltonian (constant energy shifts). 

(^^) Here we used the notation h = {h\ — hr) = 2Sz 

(22) We choose this state because it is the ground state of the Bose-Hubbard Hamiltonian in the 
absence of interactions. 
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is quantified by the number squeezing factor [159] 


(117) 




7V/4’ 


It compares the number fluctuations (AS^) to that of a coherent state (-/V/4). Therefore 
when < 1, the fluctuations are reduced compared to shot noise, and the state is said 
to be squeezed. 

One could as well define squeezing along any axes of the Bloch sphere [159]. For 
instance, phase squeezing is quantified by the phase squeezing factor = 4AS'^/iV. 

Squeezing has been shown to be a resource for interferometry below the shot noise 
limit(^^). As we saw a the beginning of section 3’2.3, the shot noise stems from the mea¬ 
surement of coherent states at the output of the interferometers. Since this corresponds 
to having no correlations between particles, or in other words, to sending the particles 
in the interferometer one after another, this limit is often considered “classical”. Inter¬ 
actions can introduce correlations between the particles, an example being squeezing, 
which can be used to obtain squeezed state at the output of the interferometer, thus 
decreasing the “shot noise” along an axis most appropriate for the measurement. 

The improvement of signal to noise ratio of an interferometer using a squeezed state 
is given, not directly by the squeezing factors introduced above, but by the so-called spin 
squeezing factor [162, 164] 


(118) 




(^ n .)2 + ( 5 „ 3 ) 2 ’ 


where = n • S and the ns are mutually orthogonal unit vectors [164]. For states close 
to the ground state of the Bose-Hubbard Hamiltonian, this expression reduces to(^^) 


(119) 


Vnas, 


It can be interpreted as the ratio of the number squeezing factor introduced earlier, 
and the coherence factor, 2{Sx)/N, characterizing the coherence of the state. Indeed, in 
the classical limit, we have 2Sx/N = cos(j), and since the mean phase is {(j>) = 0, 2{Sx)/N 
will be close to 1 when the phase distribution is peaked around 0 (very coherent state), 
and close to 0 when the phase is uniformly distributed in ] — 7r,7r] (no coherence). 


(^^) Cf. ref. [162] on the applications to spectroscopy, and references therein concerning the nse 
of squeezed states in optical interferometry. For a more general review on the nse of non-classical 
states for quantum-enhanced measurements, cf. ref. [163] and references therein. 

(^'*) This corresponds to the choice (ni,n2,n3) = (z,x,y), and noticing that that, for this 
particular state, {Sy} = 0 (or equivalently (f) = 0). We could have chosen another direction 
than 2 : to define the squeezing factor, but this choice is governed by the fact that the ground 
state is maximally squeezed along this axis. 
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Fig. 27. - Number fluctuations and coherence of the ground state of the Bose-Hubbard Hamil¬ 
tonian vs strength of interactions. The thick grey lines are obtained from an exact computation 
of the ground state of the Hamiltonian (97) with N = 100 atoms and no bias (<5 = 0). Solid 
line: dashed: coherence factor, dot-dashed: ^s- The thin solid lines are the results of the 

Bogoliubov treatment detailed in ref. [167], valid only in the Rabi and Josephson regimes. The 
dotted lines are obtained from the harmonic approximation of section 3'2.6 (see equations (130) 
and (132)). The vertical lines at 7 = 1/A^ and 'y = N separate the three regimes. 


Finally, we note that the spin squeezing factor has been related to the degree of 
entanglement present in the system [164, 165, 166]. 

3’2.4. Squeezing during splitting. We return to the splitting problem in a double-well 
potential, but now analyze how the fluctuations are modified, and show how number- 
squeezed states are generated. 

Ground-state fluctuations of the Bose-Hubbard Model. As introduced in section 3’1.4, 
double-well potentials essentially allow adiabatic splitting. The system stays in the 
ground state as long as adiabaticity is not broken by the fact that the typical time 
scale of the problem diverges when the tunnel coupling approaches J —>■ 0. It is therefore 
interesting to analyse the fluctuations of the ground state, which can be found by exact 
diagonalization of the Hamiltonian (97) as long as the particle number is not too large. 
In figure 27, the fluctuations are shown as a function of the ratio of interactions to tun¬ 
neling energy 7 = {7/2 J. The exact calculation is compared to the result of a Bogoliubov 
treatment of the Josephson junction, not detailed here [167], and to the approximate 
quadratic model detailed in section 3’2.6. 
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We see three distinct regimes: 

1. 7 <C l/N: the interactions are negligible. This is the Rabi regim, in which the 
dynamics is essentially the non-interacting dynamics presented in section 3T.3 and 
3'2.2, and the ground state is close to a coherent state. The squeezing factor is 
close to 1 (no squeezing) and the phase distribution has a typical width of l/y/N. 

2. l/N < 7 < fV: this is the Josephson regime. The dynamics was already discussed 
in section 3'2.2, in particular the small amplitude Josephson oscillations occur at 
the plasma frequency (107). In this regime the number fluctuations are reduced 
by repulsive interactions compared to the binomial case, and the coherence is also 
decreased but still high. 

3. 7 > iV: the number fluctuations are low, approaching 0 for Fock states (7 —>■ 00 ). 
This is the Fock regime. The coherence is almost completely lost. In this regime, 
we expect the mean-field treatment to be invalid because phase fluctuations are 
large. 

Squeezing during adiabatic splitting. Once again, splitting can be understood as the 
adiabatic following of the ground sate, which breaks when the adiabaticity condition 
cannot be fulfilled any more because of the divergence of 1/J{t) [168, 169]. In the many- 
body case, the Rabi frequency is replaced by the plasma frequency (107). One can show 
that the lower-lying many-body states are equally spaced by an energy fiwp in the Rabbi 
and Josephson regimes. The adiabatic condition (76) becomes: 

(120) |wpl«Wp. 

If the Josephson frequency is initially high, and splitting performed slowly enough, then 
one can reach the Josephson regime adiabatically, before adiabaticity breaks down (when 
jwpl ~ Wp). In this regime, the squeezing factors introduced above are on the order of 
iN ~ and ^0 ~ (yA^)^/"*, which corresponds to a number-squeezed state, 

because in the Josephson regime 7 > l/N (see ref. [169] and reference thereinn, and 
section3'2.6 for a derivation). 

Once adiabaticity is broken and J —> 0, no tunneling is possible and the number 
distribution is frozen. Therefore the number fluctuations cannot change anymore. The 
best spin squeezing factor attainable with this adiabatic strategy (thick grey dash-dot 
line in figure 27) has been calculated in ref. [170] to be 

( 121 ) 

which corresponds to a phase sensitivity following the Heisenberg scaling 5(j) oc 1/iV. 

Concerning the relative phase distribution, we will see in the next section that it keeps 
evolving, yielding interesting effect. 
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Fig. 28. - Evolution of the relative phase and its fluctuations in a decoupled double well, a, 
Linear evolution of the phase for various energy differences induced by tuning the angle between 
the two wells, b, Evolution of the circular standard deviation of the phase A</) corresponding 
to the orange curve of panel a. It exhibits phase diffusion. The red line is a fit to the data 
with the model of equation (123). Shaded area: theoretical prediction (see details in ref. [140]) 
c,d,e, Measured phase distributions vs time after decoupling the two wells. Figure adapted 
from ref. [140]. 


3’2.5. Phase diffusion. We now discuss the evolution of the relative phase after the 
condensate has been fully split (J = 0) in a symmetric way. In this case, the eigenstates 
of the system are the Fock states. Phase diffusion is essentially a dephasing effect between 
the Fock states involved in the superposition state describing the system [168, 159, 153, 
171]. We assume that the state after splitting is close to a coherent state (possibly 
squeezed) and therefore has small number fluctuations 2ASz = with ^ 1- The 

eigenenergy E{n) of the Fock state |(A^+n)/2, {N—n)/2) is calculated from equation (97) 
to be 

( 122 ) Ein)=E{0)-^-n+^n^. 

One can show that the linear terms leads to a global phase accumulation at the rate 
= 6t/h, and that the quadratic term, which is finite in the presence of interactions 
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Time 


Fig. 29. - Spin squeezing for optimized (symbols) and exponential splitting (gray lines). 
The quasiadiabatic exponential splitting is shown for several time scales for N=100. After 
reaches a minimum, it rises again due to loss of phase coherence caused by interaction induced 
de-phasing. With optimized splitting, better spin squeezing (lower ^r) is obtained at much 
shorter time scales. OCT calcualtions are for N=100 except: black triangle: N=500 and red 
-I- : N=1000. The horizontal lines show the Heisenberg limit ^r = Figure adapted from 

ref. [174]. 


{U 7 ^ 0), leads to a randomization of the phase. At short times, the variance of the phase 
evolves as [168, 153, 171, 169] 

(123) = 
where the phase diffusion rate is 

(124) R = 

It scales as the product of the derivative of the chemical potential with respect to the 
number of atoms ff, and the number fluctuations oc . 

Phase diffusion is not limited to Bose-Einstein condensates trapped in double-well 
potential, but rather occurs in a variety of interacting systems. It is one of the main 
factors limiting the use of dense samples for matter-wave interferometry, as the presence 
of interactions sets a limit to the coherence time of the system. Nevertheless, we stress 
that the effect can be strongly reduced by using more dilute systems, or by tuning 
interactions: either by canceling them [172, 173], or by reversing their sign in order to 
realize an analog of a spin echo in inhomogeneously broadened systems. 

A diferent approach is to optimize the dynamics of the splitting process using optimal 
control theory in order to counteract phase diffusion [175, 174, 160, 176]. The optime 
controll inspired splitting sequences of J. Grond et al. [175, 174] will allow to speed 
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up the splitting process by mre then an order of amgnitude compared to an adiabatic 
splitting and reaches the same or better number squeezing. The faster splitting leads 
to much reduced phase diffusion and consequently to much better spin squeezing ^5 
(figure 29). The simplest procedure is: first splitting fast to finite tunnel coupling, then 
keeping the double well for an extended time so that number squeezing can be establishd, 
and then splitting fast into two separated wells. 

We also underline that phase diffusion is similar to the generation of squeezing in 
the one-axis twisting scheme proposed in ref. [159] and demonstrated in the group of 
M. K. Oberthaler [177]. The variance of the phase increases, but simultaneously the 
fluctuations are reduced along other directions on the surface of the Bloch sphere, as 
squeezing is generated [159]. 

Finally, we note that at longer times, the phase will in theory display revivals [153] 
similar to those observed in optical lattices [178], and the system is expected to reach 
strongly non-classical states, for example corresponding to a superposition of two coher¬ 
ent states having a relative phase which differ by tt [161]. 

3’2.6. Effective single-particle Hamiltonian. We now show that the many-body Hamil¬ 
tonian (97) can, to some extend, be approximated by a single particle Hamiltonian. 
Similar approaches have been described in refs. [179, 170, 180]. For this we follow the 
same approach as in section 3'2.2, but keeping the population imbalance and phase as 
quantum mechanical operators(^^). From the classical limit of section 3'2.2, we have 
seen that z and (j) are canonically conjugate, therefore we are tempted to introduce the 
following commutation relation for Sz = (hr — hi )/2 and (j>: 

(125) [$,Sz]=i. 

The constant 1 x z is chosen to be dimensionless and such that the Heisenberg uncertainly 
relation yields minimal uncertainly states which resemble the coherent states introduced 
in section 3’2.3. That is ASz > 1/2, and thus = ^/N /2 ^ A^ > 1/^/N. 

We substitute I = \/hi and f = in equation (97). The Hamiltonian (97) 

becomes 

(126) H = -JN\I 1 - A{SzlNY cos li) + USl + 6Sz, 
where constant terms have again been removed. 


Harmonic approximation. Assuming no bias (5 = 0) zero imbalance {{Sz) = 0) and small 
fluctuations {ASz A, A(j) <C 1), one can develop (126) to second order in Sz and 0 to 


(^®) Note that a phase operator (f) cannot be dehned rigorously, see for instance refs. [181, 168, 
157]. 
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obtain the quadratic Hamiltonian 


(127) 


H = 




If we assume that the operators Sz and (j) have continuous spectra (which is strictly 
speaking not the case), then they are equivalent to position and momentum operators 
because they have a similar commutation relation (125), and one can therefore describe 
the many-body system by a single-particle wave function i/i(()), t). Here the phase plays 
the role of a position and the number difference, that of a momentum. In this “|()')” 
representation, the number difference operator takes the form 


(128) 


Sz = 


. d 


and the wave function satisfies the Schrodinger equation 


(129) 


■'‘a = 


2J\ 52 JN 




This is the Schrodinger equation of a single particle having a mass m (x {U + 2J/N)~^, 
in a potential having an angular frequency ojp = 2 Ji/T+lT/VT^/h. We recover the 
expression of the plasma frequency (107). 

Rabi and Josephson oscillations will therefore correspond to a dipole oscillation of the 
“phase wave function” 'tp{4>) in an effective harmonic potential. 


Ground-state fluctuations. The number and phase fluctuations can also be computed. For 
instance the ground sate is the well known Gaussian minimal uncertainty state, which 
has number and phase fluctuations given by 


(130) 

(131) 


A5, 


'^{1 + jN ) 
(l + 7jV)iA 
//V 


From the second equation, the coherence, introduced in section 3‘2.3, can be calculated 
as 


(132) 


(cos (j)) 


cos A(j) ~ 1 — 


A02 


provided that 7 <C A. The expressions (130) and (132) are represented as dotted lines in 
figure 27. We see that they agree well with both the exact calculation and the Bogoliubov 
expressions in the Rabi and Josephson regimes. 



68 


J.-F. SCHAFF, T. Langen and J. Schmiedmayer 


Phase dijfusion. Phase diffusion can also be recovered with this model. A Gaussian wave 
packet having an initial “width in momentum space” AS'^, and left in a flat potential 
(J = 0) will spread as 

(133) ^ 

which is the same expression as equation (123). 
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4. Probing Many-Body Physics by Interference 

A general understanding of quantum many-body systems is an important unsolved 
problem, touching systems from diverse fields such as cosmology, high-energy physics or 
condensed matter. In particular, the question of why and how isolated quantum systems 
relax toward equilibrium states has so far only been studied for a very limited number 
of special cases (for a review, see ref. [182]). For example, when dealing with integrable 
systems exhibiting many constants of motion, thermalization can be completely absent 
or strongly inhibited. 

A key challenge in such studies is the lack of techniques for characterizing complex 
quantum many-body states. Moreover, investigations of non-equilibrium dynamics are 
challenged by the scarcity of quantum many-body systems that are both well isolated 
from the environment and accessible to experimental study. In the following we show that 
the toolbox of atomic physics and matter-wave interferometry that has been presented 
in the previous lectures is ideally suited to tackle such problems. 

4’1. Interference of ID Bose gases. - 

4’1.1. The ID Bose gas. The experimental studies that we present are performed 
using trapped ID Bose gases. Such systems offer two unique advantages for quantum 
many-body physics: firstly, on the experimental side, realizing them with ultracold atoms 
facilitates a precise preparation and probing of the system. Secondly, on the theoretical 
side, ID Bose gases offer a model system which contains complex many-body physics but 
can still be captured with reasonable theoretical effort, in particular due to the existence 
of effective models which allow to describe the essential physics in a relatively simple 
way [183]. Furthermore, the homogeneous ID Bose gas with repulsive contact interac¬ 
tions is an example of a fully integrable quantum system [184, 185]. The approximate 
realization of such a system in experiments thus allows the study of thermalization in 
the vicinity of multiple conserved quantities and hence the study of the interplay be¬ 
tween integrability, many-body dynamics and thermalization. In this lecture we present 
a short overview of the topic. Further details can be found in the PhD thesis of Tim 
Langen [152]. 

Theoretical description. In cold atom experiments, a ID Bose gas can be realized using 
anisotropic magnetic or optical trapping potentials, where the confinement in two di¬ 
rections is strong enough such that the temperature and the chemical potential of the 
system are smaller than the excited energy levels of the trapping potential. This can be 
expressed by the condition 

(134) ksT, ij, < hujj_, 

where wj, denotes the harmonic trap frequency in the two strongly confining directions. 
In the following, we will use the convention that the strongly confining trap directions are 
the x-direction and the y-direction, such that u}± = cOx = tOy. In contrast to the strongly 
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confining directions many momentum modes can be occupied in the weakly confining 
direction. This leads to markedly different behavior than in 3D BECs, where only the 
lowest momentum mode is macroscopically occupied. These many momentum modes in 
ID Bose gases are the origin of strong density and phase fluctuations. It has been shown 
rigorously by Mermin, Wagner and Hohenberg [186, 187] that because of this enhanced 
role of fluctuations no off-diagonal long-range order and thus no BEC can exist in ideal 
ID Bose gases, even at zero temperature. 

The fact that the ID nature of the gas is achieved by strongly confining it in two 
directions of the 3D space is reflected in the scattering properties of the atoms. These 
are still 3D for the parameters reached in our experiments with ID Bose gases. For 
temperatures below the degeneracy temperature Tu = (2mkB) and sufficiently 

high density, a mean-field description is applicable to model these properties. In this 
case, the ID dynamics can be described by integrating out the two strongly confining 
directions, leading to an effective scattering potential U{z — z') = gS{z — z'). Here, 

(135) g = 2hasUj±, 


is the ID interaction parameter, with as denoting the s-wave scattering length, and we 
have assumed ^ a/ h/muj±_ [188, 189]. 

The homogeneous ID Bose gas with such delta function interactions is one of the prime 
models of mathematical and statistical physics and can be described by the Hamiltonian 


(136) 


H = 


-f - y dzdz' {z,t)'^\z',t)g5{z — z')'^{z',t)'^{z,t). 


An exact solution based on the Bethe Ansatz ]190] was found by Girardeau, Lieb and 
Liniger [184, 191, 192]. This exact solution was used by Yang and Yang to describe the 
system at finite temperature [185, 193]. Note that the existence of this exact solution 
implies that the system is integrable, which is expected to have important consequences 
on its non-equilibrium dynamics [193, 194]. In the exact solution, the strength of the 
interactions is parametrized by the Lieb-Liniger parameter 7 which is defined as 


(137) 


7 = 


mg 

K^niD 


The ID Bose gas thus becomes more strongly interacting for lower density. For 7 1 

the system is in the strongly-correlated Tonks-Girardeau regime, for 7 <C 1 it is in the 
weakly-interacting regime. The complete diagram of states at finite temperature was 
studied in ref. [195]. For the typical parameters of our experiments, 7 is on the or¬ 
der of 10“^ — 10“^ and the temperature is approximately lO^^Tp. The gas is thus 
a weakly-interacting quasi-condensate characterized by suppressed density fluctuations. 
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2 ms 
4 ms 

6 ms 

8 ms 


10 ms 


Fig. 30. - (a) Linear density profile of a ID Bose gas after a short expansion time (2 ms), 
where density fluctuations have not yet developed, (b) Linear density profile after 10 ms of free 
expansion, showing high contrast density fluctuations, where the averaged density profile (gray 
dashed line) is smooth. Studying the statistical properties of these density fluctuations enables 
the measurement of the temperature of the gas. (c) Example images of density fluctuations 
emerging during free expansion. Taken from [199]. 


The phase, on the other hand, strongly fluctuates. The corresponding correlation func¬ 
tion decays exponentially, with the thermal coherence length Xt = 2hnid/(jnksT) being 
directly related to the temperature T and the ID line density tiid of the system. 

In experiment, these phase fluctuations of a quasi-condensate can be observed in 
two different ways. First, the spatially varying phase 9{z) corresponds to a velocity 
field V = 'X/d{z) for the atoms. In time-of-flight expansion the phase fluctuations will 
thus turn into density fluctuations, similar to an optical speckle pattern [196, 197, 198]. 
An example is shown in figure 30. As the gas rapidly expands in the radial direction, 
interactions can be neglected in the expansion and it is thus possible to directly relate 
the correlation properties of these density ripples density ripples to the in situ phase 
fluctuations. Measuring these correlation properties can therefore be used to extract the 
temperature of the gas. This is particulary useful for highly-degenerate ID Bose gases, 
where no thermal background is visible anymore. Second, two quasi-condensates can 
be made to interfere, revealing a fluctuating relative phase in their interference pattern. 
This situation will be discussed in detail in the next section. 

4T.2. Multi-Mode interference. The interference of ID Bose gases in a double well can 
be described in a way that is very similar to the simple interference of the two Gaussian 
wave packets that was introduced in section 3. We also point the reader to the lectures 
of Eugene Dernier at the Varenna school 2006 [200]. To describe the many-body nature 
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of the system the fields have to replaced with operators, leading to 


= (|4'i|2) + (|4-,,|2) + 


Here, the field operators should be taken after expansion. Explicitly, ^(r, z) = t/ 2 hn 

dj 2 )-iQltl 2 hm (jgjiotes the total field operator after expansion, with Qi ^ = 
m{r ± d/2)/ht the momenta of atoms that are released from one of the condensates and 
detected at a point r. [201]. Again, the cross terms are responsible for the interference. 

To formally describe the interference pattern we introduce the operator [201, 200] 


/•^/2 

A{L)= / dz^{^ 

JL/2 


where L denotes a length scale over which the interference pattern is integrated in the 
z-direction. The complex phase of this operator can be identified with the integrated 
phase 4’{L) of the interference pattern, its magnitude |A| is connected to the interference 
contrast. For independent gases the expectation value of this operator vanishes, as the 
total phase is different in each individual realization. However, one can still observe 
high-contrast interference in these individual realizations. Consequently, the variance 


CL/2 rLj2 

(140) (|A(L)|2)=/ / dzdz' 

JLI2 Jl/2 ' 


is finite. The operator (|A(L)|^) is directly related to the mean squared contrast of the 
interference pattern {C'^{L)) = (|A(L)|^)/n^^L^, where nid = (|4'i(z)|^) is the mean 
density in the two gases. Note that the argument of the integral in equation (140) can 
be identified with the two-point correlation function of the relative phase 


C{z,z') = 


{i>rHz)'^l{z)^lHz')'^riz')) 

(|vI/,(^)P)(|4/,(zOP) ■ 


The mean squared contrast of the interference pattern is thus the double integral over 
the two-point correlation function. Neglecting the typically weak density fluctuations by 
using 'i'l r = exp(0/^r), the phase correlation function can be rewritten as 


C{z,z') = 


with (j){z) = 0i{z) — 9r{z) denoting the local relative phase. 
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Full distribution functions. Further information about the system beyond the two-point 
correlation function can be obtained from the shot-to-shot fluctuations of the contrast. 
Such noise measurements have a long and rich history in the characterization of various 
quantum systems [202]. Performing such measurements already deepened our under¬ 
standing of quantum mechanics as it led to the discovery of the Hanbury-Brown-Twiss 
effect [31, 32] which triggered the development of modern quantum optics [203]. Further¬ 
more, the study of current fluctuations led to important observations in quantum-Hall 
systems [204, 205]. Recently in atomic physics the analysis of noise correlations revealed 
the coherence properties of atom lasers [206] and enabled observations of the Hanbury- 
Brown-Twiss effect for massive fermions and bosons [207]. It was further suggested [208] 
and experimentally demonstrated [209, 210, 211, 212, 38] that noise correlations in time- 
of-flight can be used to probe strongly-correlated equilibrium states of quantum many- 
body systems. 

The shot-to-shot fluctuations of the contrast can be characterized by the moments of 
the full distribution function (FDF) P{a) of fringe contrast [213, 150] 

/|l|2m\ foo 

(143) J ' = (a™) = / P{a)a^da, 

where P{a)da measures the probability to observe a contrast in the interval between a 
and a. da. The normalized moments on the left hand side are each connected to a 
correlation function (jB]^"*) of order 2m, which is the reason why the FDF, in general, 
contains more information about the many-body state than the mean contrast. To cal¬ 
culate the full distribution function one has to calculate all moments, or equivalently, 
all even correlation functions. For the equilibrium situation powerful insights were ob¬ 
tained by mapping this problem onto a quantum impurity problem and to a generalized 
Coulomb gas model [214]. Experimentally, the FDFs have successfully been used before 
to study ID gases in thermal equilibrium [150], as well as the dynamics of an unstable 
quantum pendulum [215]. 

It is important to point out that the method of using FDFs to characterize a system 
requires the detection of single realizations of the quantum system in question. If only 
ensemble averages can be measured in the experiment, the statistics of those values will 
always be Gaussian due to the central limit theorem. Consequently, the characteristic 
higher moments of the observable will not be accessible. 

4’2. Pairs of ID Bose gases in equilibrium. - 

4’2.1. Theoretical description. A prerequisite for experiments is the ability to precisely 
prepare and characterize both non-equilibrium states, as well as the thermal equilibrium 
state of the system. One of the key advantages of coherently split ID Bose gases is that 
both these states can be controlled and described with high precision. 

For the theoretical modeling the complex many-body dynamics of the Lieb-Liniger 
Hamiltonian (equation (136) are captured using its phononic low-energy excitations. 
These excitations can be described using the Luttinger liquid formalism, which provides 
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a universal framework to describe ID quantum systems. For a single ID Bose gas one 
finds the Hamiltonian 


(144) 


Hll = 


he 

Y 


dz 



K 

TT 



The operators 6{z) and h{z) describing density and phase fluctuations are coarse-grained 
in the sense that they represent the physics in the long-wavelength limit beyond a cut¬ 
off [216]. For ID bosons the typical cutoff is defined by the inverse of the healing length 
ih = hjme. For a weakly-interacting ID Bose gas, one uses the exact solutions of the 
Lieb-Liniger Hamiltonian to obtain the analytic expressions 

(145) c=J^ , = = 

V TO y rag 

for the speed of sound c and the Luttinger parameter K. Introducing bosonic creation 
and annihilation operators Sj], and bk, the physics of the Luttinger liquid Hamiltonian 
can be rewritten as a set of uncoupled harmonic oscillators 

(146) Hll = '^ujkVJ)k, 

k^O 

where ojk = ck is the energy of the mode with momentum k. Here, we neglect the k = 0 
mode, as the spatial correlations in ID are determined by the modes with k ^ 0. The 
effect of the fc = 0 modes, on the other hand, can be identified with the phase diffusion 
that was discussed in section 3 in the context of 3D Bose-Einstein condensates. The 
excitations of the momentum modes contain both a density and a phase quadrature 

(147) bk ^ hk{z) + 9k{z). 


In the experiment we are interested in pairs of spatially separated ID Bose gases. In 
this case, the excitations contain anti-symmetric and symmetric superpositions of phase 
and density in the individual gases. These superpositions are also refereed to as the 
relative and common degrees of freedom of the system. They are given by 


(148) Yz) = 9i{z) - 9r{z) , ^com{z) = 9r{z) + 9i{z) 

for the phase, and 

(149) 9(^) = tia±M£) 


for the density. Here 9i^r{z) describes the longitudinal phase profiles and hi^r{z) the 
density of the left and the right gas, respectively. 
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The relative phase 4'{z) between the two gases determines the phase of the interference 
pattern and can therefore be measured directly in experiment. One can show that within 
the Luttinger liquid formalism the dynamics of the relative degrees of freedom and the 
dynamics of the common degrees of freedom decouple [216]. Both the dynamics of the 
relative and the common degrees of freedom is thus described by an individual Luttinger 
liquid Hamiltonian. 

Writing this Hamiltonian in momentum space and solving the Heisenberg equations 
of motion for the relative phase operator, we obtain the phase variance 


(150) (|<^/c(t)P) = sin^(a;fct) + (|(^fc(0)p) cos^(wfet). 

This calculation describes both equilibrium and non-equilibrium situations by choosing 
appropriate initial values for (|jl'fc(0)p) and (|0fc(O)p). In thermal equilibrium, the exci¬ 
tations of the system are occupied according to Boltzmann statistics = ksT/hoJk- 

In terms of density and phase quadratures this corresponds to 


(151) 

(152) 


(l^fc(O)P) 


keT 

2mkBT 
h^nidk'^' 


Assuming Gaussian fluctuations the phase correlation function in equation (142) is given 

by 


(153) C{z, z') = e-HW")--^!"')]") = g- So 


1r(l0fc(t)P>(l-COSfc|z-z'|) _ g-|z-z'|KT 


where kt = 2/Xx- As in the case of a single gas, the thermal correlation function thus 
decays exponentially, the additional factor of two being a result of the definition of the 
relative phase (j){z) in equation (148). 

4’2.2. Measurement of the matter-wave interference contrast. In a typical experimen¬ 
tal sequence we prepare a thermal ultracold Bose gas using the microscopic trap of an 
atomchip [217]. This thermal gas is then split by deforming the harmonic trap into a 
double well potential. This deformation is realized via RF dressed-state potentials [136]. 
Further evaporative cooling in the individual wells creates two completely independent 
degenerate ID Bose gases. This situation corresponds to the thermal equilibrium situ¬ 
ation. Subsequently, the gases are released from the trap and allowed to expand freely 
under gravity. After a certain time of flight they form an interference pattern, where the 
local position of the interference fringes fluctuates along the z-direction, as determined by 
4>{z). In a simple picture, every point along the length of the gas acts like an independent 
interference experiment. In every individual one of these experiments, the position of the 
fringes is determined by the local relative phase of the matter-waves emitted from the 
two wells. Stacking all these experiments together results in the observed interference 
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pattern, as shown in figure 31c. This simple picture neglects the effects of the expansion, 
but remains accurate for typical experimental parameters. 

The main experimental observables derived from this interference pattern are the 
contrast C{L) integrated over a length L and the local relative phase To extract 

them, the interference pattern is integrated along a length L. This yields a line profile 
(see figure 31) which is fitted with the cosine-modulated Gaussian function [218] 


(154) 


fL{x) = A • exp 


{x - Xo)'^\ 

2a2 J 


1 -I- C (L) cos 


/ 2Tr{x — xq) 

V 


<P{L) 


where a is the rms radius of the Gaussian profile, Xg is its center of mass, and Aj? is 
the fringe spacing. To extract the relative phase profile </)(z), the integration length L 
is set to the size of one pixel. The results are shown in figure 31. Repeated realizations 
of the experimental cycle allow us to build the time-dependent FDFs of the interference 
contrast C(L) or the phase correlation function C(z,z'). 

As in the case of a single ID Bose gas, we can further image the atoms transversally 
to extract information using the resulting density fluctuations in time-of-flight. As the 
two gases are completely independent, the resulting density ripple pattern is an inco¬ 
herent superposition of two single density ripple patterns and contains predominantly 
contributions from the symmetric mode. We simulate this situation using the Ornstein- 
Uhlenbeck process technique and use it to extract a temperature [219]. As expected for 
thermal equilibrium, we find exactly the same temperature as for the anti-symmetric 
mode. 


4’3. Probing relaxation in non-equilibrium ID Bose gases. - 
4’3.1. Goherent splitting of a ID Bose gas. 


Quenching a quantum many-body system. To study non-equilibrium dynamics we prepare 
a single degenerate ID Bose gases in a harmonic trapping potential. A quench is realized 
by rapidly deforming this harmonic trapping potential into a double well potential, as 
shown in figure 33. In general, such a quench is a rapid change of the Hamiltonian of a 
system, such that the system’s new state is no longer an eigenstate of the Hamiltonian. 
As a consequence, the system will start to evolve under the influence of the new Hamil¬ 
tonian. We are interested whether this leads to relaxation, the emergence of possible 
steady states or the emergence of thermal properties. Gurrently, such questions are un¬ 
der intense theoretical and experimental study, but in contrast to thermal equilibrium, 
no general framework exists yet to describe non-equilibrium quantum many-body sys¬ 
tems. The general problem arises from the fact that the evolution of quantum systems 
is unitary, consequently it is impossible to reach a unique thermal state from different 
non-equilibrium states. Moreover, starting from a pure non-equilibrium state, quantum 
mechanics provides no straightforward connection to mixed thermal states [182]. A pos¬ 
sible mechanism to overcome these limitations in quantum many-body systems is the 
eigenstate thermalization hypothesis [220, 221, 222]. 
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Fig. 31. - Imaging two interfering ID Bose gases, (a) After turning off all trapping potentials, 
the clouds expand and form a matter-wave interference pattern. This pattern can be imaged in 
three spatial directions. Longitudinally, the interference pattern is integrated along its length. 
In the vertical direction the local fluctuations of the interference pattern are revealed. As the 
line of sight of the vertical imaging is blocked by the atom chip, the imaging beam is reflected 
before passing the atoms. Transversally, an image of the sum of the density ripples in both gases 
can be obtained. Alternatively, the double-well trap can be turned off before the static trap, 
such that the two matter-wave packets start to roll down the potential, picking up opposite 
momenta. After time-of-flight expansion (TOF), this leads to two well-separated clouds which 
can be individually resolved using the longitudinal imaging. This procedure can thus be used 
to count the number balance of atoms in the left or right gas. (b) Emergence of the interference 
pattern in time-of-fiight in the longitudinal direction, (c) Example interference pattern (vertical 
imaging) that is used to extract phase and contrast by fitting equation (154). Figure adapted 
from [152]. 
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Fig. 32. - A pair of ID Bose gases in thermal equilibrium, (a) Two-point phase correlation 
functions C(z = z — z') for n\d = 35//im, T = (27±7) nK (red) and for n\d = 45//rm, T = (117± 
7) nK (blue). Solid lines denote the theory predictions from equation (153), including the optical 
resolution. Points are the experimental results, averaged over approximately 100 realizations. 
The temperatures used for the theory lines have been independently determined using density 
ripples, demonstrating that the gases are in equilibrium. For the hotter dataset, the correlation 
function is completely determined by the optical resolution. Such high-temperature datasets 
thus enable an independent determination of the optical resolution, (b) The corresponding 
FDFs are exponentially decaying on all length scales for hot temperatures and show a crossover 
from exponentially decaying to Gumbel-like for lower temperatures [150], both in very good 
agreement with theory [213, 219]. Figure adapted from [152]. 


In the case of the ID Bose gas the quench (i.e. the deformation of the potential) leads 
to a splitting of the gas into two uncoupled ID gases with identical phase profiles. In 
terms of excitations, this means that all the thermal excitations are initially contained in 
the common degrees of freedom. The relative degrees of freedom, on the other hand, are 
initially populated only by quantum noise created in the splitting process. The quantum 
noise adds additional energy to the system, preparing it in a non-equilibrium state. The 
aim is now to probe how the initially almost perfect correlations of the relative phase 
become obscured over time and if the thermal equilibrium state discussed above, where 
common and relative degrees of freedom are described by the same temperature, is finally 
reached [150, 224]. 
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Fig. 33. - Experimental scheme for the non-equilibrium experiments, (a) A single phase- 
fluctuating ID Bose gas is coherently split into two uncoupled parts with almost identical phase 
distributions 6i{z) and 9r{z) (phases represented by the black solid lines). These evolve for a 
variable time t. (b) At t = 0 ms, fluctuations in the local phase difference (piz) between the two 
gases are very small, but start to randomize during the evolution. The question we study is 
if, and if yes, how this randomization leads to the thermal equilibrium situation of completely 
uncorrelated gases, (c) shows typical experimental matter-wave interference patterns obtained 
by overlapping the two gases in time-of-flight (TOE). Differences in the local relative phase lead 
to a locally displaced interference pattern. Integrated over a length L, the contrast C{L) in 
these interference patterns is a direct measure of the strength of the relative phase fluctuations 
and thus enables the investigation of the dynamics. Figure taken from ref. [223]. 


4’3.2. Dynamics of the matter-wave interference contrast. To this end, the two gases 
are allowed to evolve in the double-well potential for a varying evolution time t before 
the relative phase correlations are probed via time-of-flight matter-wave interference 
(Fig 33c). As described above, increasing fluctuations in the relative phase lead to a 
locally displaced interference pattern, making matter-wave interferometry a direct probe 
for the dynamics of the system [201, 225, 214, 137, 226]. 

In figure 34a we show the result of this procedure for different integration lengths and 
different evolution times after the splitting [218]. We observe a rapid, length-dependent 
decay of the contrast on time-scales of 5 ms {L — 18 /rm integration length) to 20 ms (L = 
100 fim integration length). These time-scales are much faster than the ones expected for 
thermalization, which is thought to be very slow in ID Bose gases [227]. The observations 
demonstrate that the system relaxes but they do not reveal the physical mechanism 
behind the relaxation. To reveal the nature of the quasi-steady state, we go beyond mean 
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Fig. 34. - Dynamics of a coherently split ID Bose gas. (a) Measured values of the mean squared 
contrast for various integration lengths (points). From top to bottom: L = 18,40,60,100^m. 
The lines show the results of a Luttinger liquid calculation for these integration lengths [218, 
216]. We observe a relaxation process in which a steady-state is established on a time-scale 
depending on L. (b) Experimental non-equilibrium distributions (histograms) of the matter- 
wave interference contrast for this steady state. The solid red lines show theoretical equilibrium 
distributions with an effective temperature of Tefi = 14 nK, which is significantly lower than 
the true initial temperature of the gas (T = 120nK). The prethermalized nature of the state 
is clearly revealed by comparing it to the vastly different thermal equilibrium situation shown 
in (c), which can be prepared by creating two completely independent ID Bose gases. Figure 
adapted from [223]. 


values and measure the full distribution function (FDF) of the contrast. FDFs describing 
the evolution are shown in figure 35. They reveal distinctively different behavior on 
different length scales, again demonstrating the multimode nature of ID Bose gases. 
The results for the steady state are shown in figure 34b. For long integration lengths the 
distributions decay exponentially, reflecting the randomization of the phase and the decay 
of coherence. On short integration lengths, however, the histograms show a Gumbel-like 
distribution with a large probability for high contrasts. This reveals that some coherence 
from the initial state still remains in the system. The quasi-steady state can directly be 
compared with the thermal equilibrium state, which can be prepared experimentally by 
splitting a thermal gas and cooling the resulting two gases independently to quantum 
degeneracy. The resulting FDFs are shown in figure 34c. They decay exponentially for 
all length scales, as expected for a thermal equilibrium state in our temperature range, 
which typically exhibits very little coherence. From this strong difference, we conclude 
that the dynamically emerging quasi-steady state is not the thermal equilibrium state of 
the system. 

4’3.3. Prethermalization. Although it is not the thermal equilibrium state, the quasi¬ 
steady state still shows thermal properties. This is revealed by fitting the experimental 
data with the theoretically expected shape for a thermal equilibrium state of a certain 
temperature. For the experimentally prepared thermal equilibrium state this procedure 
yields, as expected, a temperature which corresponds to the temperature determined 
from the density ripples. In contrast to this, the quasi-steady state reveals thermal 
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Fig. 35. - Multimode dynamics revealed by PDFs of matter-wave interference. The probability 
density of contrast C{L) and phase ^{L) of interference patterns is measured for each integration 
length L (horizontal axis) and evolution time t (vertical axis). For each realization, they are 
plotted in a polar plot, where the radius denotes the contrast and the angle corresponds to the 
phase. Averaging over many realizations results in a probability distribution, where red and 
blue denote high and low probability respectively. For each value of L, the right (left) columns 
correspond to the experimental data (theoretical calculations), and full cloud to L = 100 fim. At 
t = 1.5 ms the high contrasts and small phase spreads demonstrate the coherence of the splitting 
process. As time evolves, a steady state emerges and two distinct length-scale dependent regimes 
appear: a regime where only the phase diffuses (A) and a regime where the contrast decays 
strongly (B). For short (long) L, the phase is random and the probability of observing a high 
contrast is high (low), resulting in a ring (disk) shape in the density plot. Figure from ref. [218]. 


PDFs corresponding to an effective temperature that can be significantly lower than the 
temperature determined from the density ripples. For the data presented in figure 34b we 
find Teff = 14 nK, almost an order of magnitude lower than the temperature T = 120 nK 
of the initial gas. 

Emergence of an effective temperature. The nature of this very low effective temperature 
can be understood in the following way. The PDFs measure the temperature of the 
anti-symmetric (or relative) degrees of freedom. As outlined above, these degrees of 
freedom are created during the splitting process and only contain quantum noise. This 
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quantum noise that is introduced into the relative degrees of freedom has its origin in the 
binomial distribution of atoms into the left and right condensates during the splitting. 
We assume this splitting process to be instantaneous, i.e. tgpiit ^ fi//i = (Jh/c, so that 
no information can propagate along the length of the system. Therefore the splitting of 
the individual atoms is completely uncorrelated, resulting in binomial statistics of the 
splitting process [216, 228]. For a segment with a length (Jh containing N = 2nid<Jh 
atoms, this leads to relative atom number fluctuations that are given by 7V/4. In terms 
of the density quadrature this can be expressed as 

(155) (l^>fc(0)P) = 

The phase quadrature follows from the uncertainty relation = —i5k,k' as 

(156) {\M0)\^) = 

2nid 

For the excitation numbers this corresponds to 

(157) (515,) = 

where ksTeS = gnid/^. is the energy that is added to the system in the splitting pro¬ 
cess. These occupation numbers are vastly different than their counterparts in thermal 
equilibrium (see equation (152)). In detail density fluctuations are strongly enhanced, 
whereas phase fluctuations are strongly suppressed. However, as in thermal equilibrium, 
all momentum modes are initialized with the same amount of energy, independent of 
their momentum. Dephasing of the many momentum modes due to their different fre¬ 
quencies ujk then leads to the thermal-like statistics of the FDFs. Probing the system on 
different length scales acts as a filter for the effect of the different modes. For a given 
integration length L only modes with a wavelength shorter than L can contribute to the 
randomization of the phase. Thus the level of contrast in the quasi-steady state is high 
for short L, but significantly lower for long L, where many modes can contribute. 

Evidence for a generalized Gibbs state. The symmetric degrees of freedom of the system 
are still hot, because they are thermally populated from the beginning. This observa¬ 
tion has two important consequences. Firstly, anti-symmetric and symmetric degrees of 
freedom have not equilibrated and the system has thus not thermalized. Secondly, the 
quasi-steady state can be identified with a prethermalized state, as predicted for heavy- 
ion collisions [229], condensed matter systems [230, 231] or split ID Bose gases [216]. 
In contrast to usual statistical mechanics, more than one temperature is needed to de¬ 
scribe the system. This suggests that the relaxed state is an example of a generalized 
Gibbs state [232], which is generally believed to describe relaxed states in systems with 
dynamical constraints, such as our ID Bose gas. 
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4’3.4. Light-cone-like emergence of thermal correlations. More insights insights into 
the emergence of the prethermalized state can be obtained via time-resolved measure¬ 
ments of the two-point phase correlation function. The results of these measurements 
are shown in figure 36. Directly after the quench, the phase correlation function C{z,t) 
is close to unity for any distance z = z — z', where z and z' are arbitrary points along 
the system. This is a direct manifestation of the long-range phase coherence produced 
by the splitting process.After a given evolution time t, the phase correlation function 
decays exponentially up to a characteristic distance Zc and stays nearly constant after¬ 
wards: C{z > Zc, t) = C{zc, t). This means that beyond the distance Zc long-range phase 
coherence is retained across the system. With longer evolution time, the position of Zc 
shifts to larger distances and the value oiC{z > Zc, t) gradually decreases. The evolution 
continues until the system reaches the prethermalized state, where the correlations decay 
thermal-like throughout the entire system as 

(158) C(z,z') = 

where Xxeff = fi^nid/rnkBTeff is the prethermalized correlation length [218]. 

Extracting the crossover points Zc we observe a clear linear scaling of the position 
Zc = 2ct, characterizing the local decay of correlations with time. This observation re¬ 
veals that an arbitrary point in the gas loses its correlations with other points up to 
a certain separation Zc, while long-range phase coherence persists outside this horizon. 
The experimental data thus show that the prethermalized state locally emerges in a light- 
cone-like evolution, where c plays the role of a characteristic velocity for the propagation 
of correlations in the quantum many-body system. For the data presented in hgure 36b 
a linear fit allows to extract a velocity of c = 1.2 ± O.lmm/s, which is in good agree¬ 
ment with a Luttinger liquid calculation including the trapping potential [233]. This 
calculation reveals that the underlying mechanism is the spreading of excitations in the 
system. These excitations can be pictured as quasi-particles which transport information 
about the quench through the system with their caracteristic velocity c [234, 235, 236]. 
If the quasi-particles emitted from two points separated by Zc meet after t = zd^c these 
points establish thermal correlations. This provides evidence for the local relaxation 
conjecture [237], which provides a general link for the spreading of correlations and the 
establishment of thermal properties in quantum many-body systems. 

4’4. Conclusion. - Matter-wave interferometry with ID Bose gases offers important 
insights into the world of quantum many-body systems. For example, our examples 
demonstrate that non-equilibrium dynamics are far richer than a simple direct relaxation 
to thermal equilibrium. The pairs of ID Bose gases introduced in this lecture also allow 
the realization of many more complex non-equilibrium situations to be explored, for 
example by introducing a tunnel coupling between the two gases. We expect matter- 
wave interferometry therefore to have profound implications for our understanding of 
the emergence of thermal and classical properties in quantum systems, the study of 
which is an ongoing theoretical and experimental endeavor. 
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Fig. 36. - Local emergence of thermal correlations in a light-cone-like evolution, (a) Experi¬ 
mental phase correlation functions C(z,t) (filled circles) compared to theoretical calculations 
(solid lines). From top to bottom, the evolution time t increases from 1ms to 9ms in steps 
of 1 ms. The bottom (green) line is the theoretical correlation function of the prethermalized 
state. For each t, the constant values of C{z, t) at large « can be used to determine the crossover 
distance Zc{t) up to which the system forgets the initial long-range phase coherence (see text 
for details), (b) Position of the crossover distance Zc as a function of evolution time t, revealing 
the light-cone-like decay of correlations. Error bars denote the uncertainty in Zc, following from 
the standard deviation of the constant values of C{z, t) and the uncertainty in the effective tem¬ 
perature of the prethermalized state (see Supplementary Information). The solid line is a linear 
fit, the slope of which corresponds to twice the characteristic velocity of correlations. Inset: 
schematic visualization of the dynamics. The decay of correlations is characterized by a front 
moving with a finite velocity: for a given time t, C(z,t) is exponential (thermal) only up to the 
characteristic distance Zc{t). Beyond this horizon, long-range phase coherence is retained. Note 
that in the experimental data shown in (a), the sharp transitions are smeared out by the finite 
experimental imaging resolution. Eigure adapted from [238] 
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